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ABSTRACT 


The  total  energy  system  concept  has  been  proposed  as  a 
possible  means  of  reducing  the  cost  of  providing  electricity 
at  MIT.   An  overview  of  key  factors  influencing  the  possible 
shift  to  a  total  energy  system  approach  is  presented.   Campus 
steam  and  electrical  load  profiles  are  defined  and  the  depen- 
dence of  load  upon  ambient  temperature  is  analyzed.   Load 
growth  and  the  future  impact  of  conservation  measures  at  MIT 
are  addressed  in  relation  to  the  relative  sizing  of  a  proposed 
total  energy  plant.   A  demand  model  is  constructed  for  use  in 
simulating  the  operation  of  alternative  total  energy  designs 
on  a  computer.   A  comparison  of  1976  consumption  data  at  MIT 
with  that  predicted  by  the  load  model  is  made,  establishing 
the  validity  of  the  model  for  further  use  in  total  energy 
system  simulation.   Methods  of  modeling  different  equipment 
configurations  are  discussed  for  the  purpose  of  devising 
computer  programs  to  aid  in  comparative  cost  studies. 
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I    INTRODUCTION 

The  cost  of  providing  for  MIT's  energy  needs  [1] 
increased  from  $1.8  million  in  fiscal  year  1970  to  $5.1 
million  in  fiscal  year  1974.   The  last  three  years  alone  have 
seen  a  73%  increase  in  the  annual  cost  of  electricity  which 
is  purchased  through  the  local  utility  (Cambridge  Electric) . 
Viewed  by  themselves,  these  figures  hardly  seem  startling  as 
they  reflect,  at  the  very  least,  the  sharp  rise  in  fuel 
prices  attendant  to  the  1973  oil  embargo.   Cost  figures  alone, 
however,  distort  the  picture  of  energy  consumption  at  MIT, 
for  these  monetary  increases  have  occurred  in  spite  of  signi- 
ficant energy  conservation  efforts.   More  specifically,  in 

the  past  three  years  the  intensity  of  electricity  consumption 

2 

at  MIT  (measured  in  KWH/ft  -year)  has  dropped  by  23.5%. 

Although  energy  conservation  measures  continue  at  MIT  and, 
in  all  likelihood,  will  further  reduce  the  average  kilowatt 
load,  the  cost  of  providing  electricity  is  certain  to  keep 
increasing.   These  facts  provide  the  motivation  behind  a  study 
to  determine  more  cost  effective  means  of  supplying  energy 
to  the  MIT  campus. 

Implicit  in  the  above  is  the  concept  of  on-site  genera- 
tion of  electricity.   To  what  extent  this  might  involve 
divorcing  MIT  from  its  present  utility  ties  is  a  question 
which  ultimately  could  determine  the  feasibility  of  such  an 
undertaking.   Nonetheless,  if  annual  costs  are  to  be  moderated, 
some  level  of  electrical  generation  is  needed.   Recognizing 
that  MIT  has  sizable  thermal  loads  the  year  around  (both 
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heating  and  air  conditioning) ,  the  on-site  generation  of 
electricity  translates  to  "total  energy  system". 

Offering  the  advantage  of  operating  efficiencies  in  the 
range  of  65-80%,  total  energy  systems  make  efficient  utiliza- 
tion of  the  thermal  energy  which  is  a  necessary  by-product  of 
electrical  power  generation  [2],   It  is  conceivable  that  with 
the  installation  of  electrical  generators  and  their  associated 
prime  movers,  all  of  MIT's  power  needs  could  be  provided  from 
one  fuel  source.   Typically,  one  third  of  the  available  fuel 
energy  is  used  to  produce  electricity  in  any  power  generation 
scheme.   The  more  efficient  the  recovery  of  the  remaining  two 
thirds  available  energy,  the  more  economically  justifiable  is 
the  chosen  total  energy  design.   As  total  energy  systems  tend 
to  exhibit  higher  first  costs,  their  attractiveness  lies 
solely  in  the  lower  annual  operating  costs  they  can  provide. 
On  the  surface,  therefore,  total  energy  system  schemes 
warrant  investigation  to  determine  their  cost  relative  to 
the  methods  presently  employed  to  provide  MIT's  energy  needs. 

A  study  such  as  this  must  necessarily  begin  with  a 
thorough  assessment  of  thermal  and  electrical  loads  at  MIT. 
The  more  accurate  this  evaluation,  the  more  tailored  the 
specific  total  energy  design  will  be  to  meet  the  required 
demands.   Unlike  engineering  consulting  firms  which  normally 
must  estimate  loads  in  new  buildings  for  the  purposes  of  power 
plant  sizing,  MIT  is  fortunate  to  have  available  detailed 
steam  and  electrical  data  from  which  load  profiles  may  be 
constructed.   An  unnecessary  source  of  error  is  thereby 
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eliminated.   Considerable  effort  has  been  devoted  to  this 
task  with  the  result  that  comprehensive  daily  electrical 
and  steam  usage  profiles  are  available  for  all  seasons 
(Chapters  III  &  IV) . 

Two  areas  which  impact  heavily  on  the  analysis  of 
Institute  loads  are  long  term  building  plans  and  the  future 
direction  of  campus  conservation  measures.   They  are 
addressed  in  Chapter  V.   For  the  purposes  of  system  design 
and  selection,  it  has  been  assumed  that  any  power  plant  con- 
figuration must  accommodate  projected  load  growth  to  the 
year  1990.   This  requirement  insures  against  gross  under- 
specification  of  plant  capacity  to  meet  those  demands.   By 
quantifying  the  effects  of  future  campus  energy  conservation 
(above  and  beyond  those  measures  which  have  already  been 
taken) ,  the  overall  design  capacity  of  a  proposed  plant  may 
be  scaled  down  somewhat  over  what  it  might  otherwise  be. 

The  conservation  measures  of  interest  are  those  which 
constitute  the  Facilities  Management  System  [3].   Implemented 
officially  in  the  fall  of  1976,  this  program  is  concerned  with 
conservation  through  power  management.   Under  the  control  of 
dual  PDP-11/40  processors,  virtually  every  building  on 
campus  will  eventually  have  its  HVAC  systems  regulated  by 
preset  on/off  commands.   This  will  serve  to  decrease  power 
consumption  by  automated  equipment  shutdown  during  periods 
when  building  usage  does  not  justify  full  scale  equipment 
operation.   Taking  into  account  this  reduction  of  campus 
loads  afforded  through  the  FMS  could  obviate  the  need  for 
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designing  a  plant  to  accommodate  full  scale  electrical  HVAC 
usage  at  night. 

The  load  profiles  which  have  been  obtained,  in  conjunc- 
tion with  information  pertinent  to  campus  growth,  permit  a 
determination  of  an  upper  bound  for  required  plant  capacity. 
Decisions  regarding  what  portion  of  the  electrical  load  MIT 
might  assume  will  be  dictated  by  the  results  of  cost  trade- 
off studies  with  Cambridge  Electric  to  determine  acceptable 
rate  structures  for  the  purchase  of  supplementary  power  to 
the  campus.   As  a  means  of  facilitating  this  exchange,  a 
methodology  has  been  outlined  to  permit  the  computer  modeling 
of  certain  total  energy  system  design  schemes.   For  any  level 
of  power  generation  required,  the  specific  plant  equipment 
capacity  may  be  easily  modified  so  as  to  accommodate,  in  the 
most  efficient  manner,  the  requisite  thermal  loads.   The  main 
computer  program,  which  is  described  more  fully  in  Chapter  VI, 
uses  as  input  data  load  information  which  is  representative 
of  a  typical  year  at  MIT.   It  passes  daily  kilowatt  and  steam 
demands  to  separate  numerical  simulation  subroutines,  each  of 
which  models  a  specific  total  energy  system  alternative  (gas 
turbine,  steam  turbine,  diesel) .   The  user  is  provided  with 
output  in  the  form  of  annual  operating  cost  information  for 
the  chosen  equipment  configuration. 

The  problem  MIT  faces  regarding  the  possible  shift  to 
total  energy  is  a  complex  one.   It  involves  far  more  than  the 
mere  selection  of  equipment  for  a  new  power  plant.   Underlying 
any  decision  about  capital  expenditure  is  the  realization  that  a 


-17- 

sizable  power  plant  now  already  exists  -  one  which  cannot 
simply  be  scrapped  but  which  in  some  manner  must  be  integrated 
into  a  more  efficient  arrangement  for  power  generation. 
Additionally,  subjects  such  as  fuel  availability  and  environ- 
mental restrictions  must  be  addressed  in  their  entirety  as 
they  could  exert  strong  influence  on  the  decision  making 
process.   It  is  not  the  intention  herein  to  investigate  in 
sufficient  detail  all  the  outlying  factors  which  must  be 
considered  prior  to  final  plant  selection.   Indeed,  the  bulk 
of  what  follows  is  concerned  only  with  load  modeling.   By  way 
of  placing  this  thesis  in  proper  perspective,  however,  Chapter 
II  has  been  included.   It  summarizes  pertinent  information 
on  the  existing  plant  facility  and  addresses,  in  brief  form, 
areas  for  later  investigation  as  part  of  the  overall  total 
energy  study. 
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II   SUPPLEMENTARY  INFORMATION 
2. 1   Background 

Decisions  relevant  to  the  selection  of  "candidate" 
total  energy  systems  for  MIT  must  be  made  with  a  thorough 
knowledge  of  existing  facilities.   It  is  appropriate,  there- 
fore, to  review  the  present  arrangement  both  for  the  genera- 
tion of  steam  and  the  purchase  of  electricity  from  Cambridge 
Electric. 

2.1.1   Steam  Generation 

The  Central  Utility  Plant  houses  five  boilers. 
Four  of  these  units  stand  in  the  original  building,  con- 
structed in  1916.   The  fifth  unit  lies  to  the  west  of  the 
other  four  in  a  building  extension  which  has  enough  space  to 
allow  for  doubling  the  capacity  of  the  units  now  in  the 
original  building  [4].   Altogether,  the  installed  capacity 
of  the  boilers  is  400,000  lbs/hour.   Steam  is  generated  at 
200  psig,  425°F.   All  turbine  driven  auxiliary  equipment  uses 
steam  at  these  conditions.   In  addition,  the  four  steam  tur- 
bine driven  centrifugal  compressors  for  the  chiller  plant 
use  this  steam.   Total  capacity  of  this  plant  is  10,500  tons. 
These  four  steam  turbines  are  straight  condensing  and  utilize 
cooling  towers  to  effect  the  heat  transfer  necessary  to  ensure 
a  steady  supply  of  low  temperature  water  for  their  operation. 

Exhaust  steam  from  the  turbine  driven  auxiliaries  is 
provided  at  a  pressure  of  5  psig.   A  common  20  inch  header 
distributes  this  steam  to  one  portion  of  the  Institute's 
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main  building  group  for  hot  water  and  space  heating  purposes. 
When  the  heating  demand  for  this  area  cannot  be  met  by  normal 
exhaust  steam,  augmenting  steam  is  provided  to  the  20  inch 
header  through  a  200-5  psig  reducing  station.   For  the 
remainder  of  the  campus,  200  psig  steam  is  distributed  directly 
to  individual  buildings  where  it  is  reduced  in  pressure 
locally  for  heating  purposes. 

The  Central  Utility  Plant  is  designed  to  utilize  either 
low  sulfur  content  oil  or  natural  gas  as  its  fuel.   In  the 
New  England  area  the  primary  source  of  energy  for  industrial 
users  has  in  recent  years  been  #6  residual  fuel  oil.   Most  of 
this  is  imported.   Fuel  oil  storage  is  divided  among  three 
locations  on  campus.   The  combined  storage  capacity  is  pre- 
sently 550,000  gallons.   This  corresponds  to  a  winter  reserve 
under  current  loads  of  approximately  two  to  two-and-one-half 
weeks. 

Modernized  to  include  electronic  automatic  combustion 
controls,  the  present  Central  Utility  Plant  has  an  average 
boiler  operating  efficiency  of  approximately  83%.   A  summary 
of  the  now  existing  facility  is  provided  below. 


INSTALLATION 

EQUIPMENT 

DESIGN  CAPACITY 

70,000  lb/hr. 
70,000  lb/hr. 
80,000  lb/hr. 
80,000  lb/hr. 
100,000  lb/hr. 

DATE 

Boiler  No.  1 
2 
3 
4 
5 

1950 
1950 
1964 
1964 
1971 

'ig.  Unit  No. 

1 
2 
3 
4 

1,500  tons 
1,500  tons 
3,500  tons 
4,000  tons 

1967 
1967 
1973 
1975 

•20 


EQUIPMENT 

Auxiliary  Machinery 

Feed  Pumps 

2  motor  driven 
2  turbine  driven 


Fuel  Oil  Pumps 
1  motor  driven 
1  turbine  driven 


DESIGN  CAPACITY 


18,070  gal/hr 
18,070  gal/hr 


3,240  gal/hr 
3,240  gal/hr 


INSTALLATION 

DATE 

1965  -  1974 

1965  -  1974 

1965 

1965 

2.1.2   Electricity  Supply 

Since  1938,  when  MIT  first  entered  into  a 
purchase  agreement  with  Cambridge  Electric  Light  Company, 
the  Institute  has  relied  upon  central  station  power  generation 
almost  exclusively.   In  1972,  a  925  KVA  diesel  engine  genera- 
tor was  added  to  the  Central  Utility  Plant  for  the  purpose  of 
peak  shaving.   It  presently  operates  on  a  somewhat  irregular 
basis,  serving  to  reduce  the  billing  peaks  during  the  hours 
0800  -  1800  on  weekdays. 

Electricity  from  Cambridge  Electric  is  fed  to  the 
Institute  through  a  system  of  primary  loops.   At  present  this 
system  consists  of  three  13.8  KV  switching  stations,  each 
having  two  incoming  feeders.   One  of  the  feeders  is  common 
to  two  of  the  switching  stations.   Each  switching  station 
distributes  power  to  the  campus  through  an  arrangement  of 
looped-primary  feeders.   The  system  serves  a  major  part  of 
the  campus  loads  directly  at  13.8  KV.   It  also  feeds  four 
transformers  which  supply  a  2300  V  distribution  system. 

MIT  is  presently  billed  for  electicity  under  the  Rate-8 
structure.   Demand  charges  are  based  on  the  peak  metered 
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kilowatt  load  during  each  30  day  billing  period.   A  15.6% 
rate  adjustment  is  presently  in  effect  under  Rate-8.   In 
addition,  a  fuel  adjustment  charge  of  2.549C  per  kilowatt 
hour  has  been  imposed.   A  recent  cost  figure  for  purchased 
electricity  at  MIT  is  3.6C  per  kilowatt  hour. 

2. 2   Problem  Overview 

The  importance  of  this  particular  study  to  the  task 
of  determining  more  cost  effective  means  of  supplying  campus 
energy  cannot  be  overemphasized.   It  is  a  necessary  and  pur- 
posefully comprehensive  "first  step"  toward  the  possible 
adoption  of  a  scheme  for  self-generation  of  electrical  power 
at  MIT.   It  should  not,  however,  be  viewed  as  all  encompass- 
ing.  Even  after  the  results  of  this  analysis  are  presented, 
serious  questions  will  remain  concerning  the  implementation 
of  any  total  energy  system  at  MIT. 

2.2.1  Arrangements  With  Local  Utility 

A  foreseeable  trend  is  developing  within  the 
utility  sector  of  the  United  States  -  one  which  seems  certain 
to  win  approval  of  the  present  Administration.   This  trend  is 
toward  the  peak  load  pricing  of  electricity.   Very  simply, 
present  utility  pricing  policy  discourages  conservation. 

Incorporated  in  the  President's  National  Energy  Plan, 
which  was  submitted  to  Congress  April  29,  1977,  are  recommen- 
dations for  sweeping  utility  reform  legislation  [5].   It  has 
been  proposed  that  electric  utilities  be  required  to  offer 
daily  off-peak  rates  to  each  customer  who  is  willing  to  pay 
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metering  costs  or  provide  a  direct  load  management  system. 
MIT  essentially  fits  into  both  these  categories.   More 
important,  however,  are  the  statements  make  by  the  President 
about  cogeneration. 

The  simultaneous  production  of  process  steam  and 
electricity,  cogeneration  is  simply  another  word  for  total 
energy.   At  present,  a  variety  of  institutional  barriers 
impede  its  development  for  wide  scale  use  in  industry.   Chief 
among  these  is  an  almost  uniform  resistance  on  the  part  of 
local  utilities  to  allow  their  lines  to  run  in  parallel  with 
total  energy  lines.   This  could  prevent  the  application  of 
cogeneration  in  building  complexes  which  have  something  other 
than  an  equitable  mix  of  thermal  and  electric  loads  -  situa- 
tions where  a  utility  tie-in  could  be  economically  advanta- 
geous.  Another  barrier  to  development  of  cogeneration  schemes 
is  the  comparatively  high  first  cost  to  the  builder.   As  this 
necessitates  long  term  investment  in  order  to  satisfy  the 
more  lengthy  amortization  periods  for  financing,  strong 
investment  incentives  must  exist  before  a  major  committment 
of  capital  funds  is  made.   In  the  past  these  incentives  have 
been  limited  to  the  ultimate  life-cycle  cost  savings  associ- 
ated with  on-site  generation  of  electricity  as  opposed  to 
purchase  from  a  utility.   In  view  of  the  risks  involved  in 
such  an  undertaking,  however,  it  is  certain  that  government 
must  make  cogeneration  more  attractive. 

Citing  1975  statistics  which  show  that  waste  heat  in  the 
industrial  and  utility  sectors  accounts  for  over  7  million 
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barrels  of  oil  per  day  in  the  U.S.,  President  Carter  has 
outlined  a  rather  comprehensive  program  to  encourage  cogenera- 
tion  in  the  National  Energy  Plan  [5],   It  has  been  proposed 
that  firms  generating  their  own  electricity  be  assured  of 
receiving  fair  rates  from  utilities  for  both  the  surplus 
power  they  might  sell  and  for  the  backup  power  they  might  buy. 
Moreover,  the  President  has  suggested  that  industries  using 
cogeneration  (MIT  would  fit  into  this  category)  be  exempt 
from  State  and  Federal  public  utility  regulation.   In  addi- 
tion, they  would  be  entitled  to  use  public  utility  trans- 
mission facilities  to  sell  surplus  and  purchase  backup  power. 
By  way  of  easing  the  first  cost  to  the  builder,  an  additional 
tax  credit  of  10%  above  the  existing  investment  tax  credit 
is  proposed  for  cogeneration  equipment.   A  key  feature  of 
the  Energy  Plan,  and  one  which  could  greatly  facilitate  the 
selection  of  a  total  energy  system  for  MIT  is  the  provision 
whereby  industrial  firms  which  invest  in  cogeneration  equip- 
ment could  be  exempt  from  the  requirement  to  convert  from 
oil  and  gas  in  cases  where  the  exemption  is  necessary  for 
cogeneration.   This  is  particularly  significant  in  the  New 
England  area  where  coal  is  presently  not  available  commer- 
cially in  sufficient  quantities  to  sustain  a  widespread 
application  of  coal-based  total  energy  systems. 

It  seems  clear  that  the  future  of  total  energy  in  this 
country  is  bright.   The  present  Administration  is  rapidly 
paving  the  way  for  increased  acceptance  by  local  utilities  of 
cogeneration  schemes.   What  form  the  government  regulations 
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will  ultimately  take  is  open  to  speculation.   Indeed,  final 
resolution  of  the  many  questions  surrounding  utility  rate 
reform  is  possibly  years  away.   It  behooves  MIT,  therefore, 
to  follow  closely  the  developments  in  this  area.   Decisions 
pertinent  to  total  energy  system  selection  must  be  constantly 
reviewed  and,  if  needed,  revised  in  response  to  the  perceived 
changes  in  utility  rate  structure,  provisions  for  purchase 
agreements,  and  governmental  incentives  for  equipment  invest- 
ment. 

2.2.2   Waste  Heat  Mangement 

For  any  total  energy  system  to  be  successful 
the  facilities  it  services  must  have  a  reasonably  steady 
thermal  demand  in  relation  to  the  power  generated.   This  is 
partially  the  case  with  MIT  as  it  is  characterized  by  a 
strong  heating  load  during  winter  months  and  employs  on  a 
regular  basis  steam  driven  air  conditioning  compressors  during 
the  summer.   While  the  magnitudes  of  electrical  demand  are 
not  vastly  different  from  one  season  to  another,  those  of 
steam  demand  are. 

The  average  winter  daily  steam  load  is  approximately 
twice  that  of  the  summer.   This  suggests  that  whatever  total 
energy  schemes  are  proposed  should  be  leveled  at  satisfying 
the  requisite  summer  electrical  and  thermal  loads  under 
normal  operation.   The  excess  heat  required  in  the  winter 
most  probably  will  come  from  augmenting  the  generation  of 
steam  in  some  fashion.   The  management  of  waste  heat  on  a 
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daily  basis,  however,  is  considerably  more  involved  than 
this. 

Waste  heat  management  concerns  the  storage  of  heat  which, 
because  of  an  imbalance  in  hourly  thermal  and  electric  loads, 
exists  as  a  by-product  of  electrical  power  generation.   Once 
the  prime  movers  of  a  total  energy  system  are  chosen,  detailed 
calculations  must  be  performed  to  determine  what  heat  recovery 
system  should  be  employed.   A  lengthy  analysis  in  itself, 
the  proper  resolution  of  waste  heat  allocation  impacts  greatly 
on  system  feasibility  studies.   It  requires  a  thorough  assess- 
ment of  the  magnitude,  duration  and  coincidence  of  electrical 
and  thermal  loads  [ 6]  for  the  purpose  of  determining  the 
"worst  possible  mix"  of  the  two. 

Although  not  addressed  specifically  in  this  study,  the 
sizing  of  waste  heat  storage  devices  is  an  integral  part  of 
the  preliminary  design  and  analysis  of  a  total  energy  system. 
It  necessarily  must  follow  the  initial  equipment  selection 
phase  and  must  provide  feedback  information  on  possible  alter- 
native equipment  choices. 

2.2.3   Role  of  Renewable  Resources 

It  is  highly  unlikely  that  either  wind,  wave 
or  geothermal  power  will  ever  be  used  to  supplement  the  energy 
needs  at  MIT.   Neither  the  Institute's  size  nor  location  will 
permit  it.   From  the  standpoint  of  future  energy  needs, 
though,  there  is  a  distinct  possibility  that  solar  insolation 
will  play  a  role. 
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Based  on  research  which  has  been  conducted  at  the 
University  of  Delaware  [2],  mass  production  of  reasonably 
efficient  thin  film  photovoltaic  cells  could  become  competi- 
tive with  central  station  power  generation  by  the  end  of  the 
1980 's.   Although  this  is  a  matter  of  considerable  debate 
within  the  solar  energy  community,  there  is  little  doubt  that 
such  proposals  will  receive  the  increasing  attention  of 
policy  planners  in  years  to  come.   Currently,  photovoltaic 
systems  are  economic  only  for  small  decentralized  applications; 
however,  the  potential  for  price  reductions  which  would  make 
them  economical  for  a  broader  range  of  applications  is 
dramatic. 

The  most  likely  manner  in  which  solar  energy  will  be 
used  at  MIT  is  in  the  heating  and  cooling  systems  of  new 
buildings,  perhaps  dormitories.   Such  usage  has  been  demon- 
strated feasible  in  other  areas  of  the  country,  notably 
Texas,  where  an  entire  entension  of  the  North  Lake  Community 
College  [7]  has  been  designed  as  a  solar  total  energy  system. 
More  convincing  evidence  that  savings  can  be  achieved,  however, 
is  needed  for  the  New  England  area.   In  that  sun  cover  in 
Boston  averages  only  55%  over  a  year,  several  prototype 
installations  are  required  in  this  region  to  determine  system 
feasibility. 

The  results  of  solar  demonstration  programs  now  being 
carried  out  by  the  Energy  Research  and  Development  Admini- 
stration and  the  Department  of  Housing  and  Urban  Development 
[  5  ]  will  help  provide  some  of  the  much  needed  information 
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about  solar  product  reliability.   Similarly,  the  proposed 
installation  of  solar  equipment  in  federal  office  buildings 
will  serve  as  a  basis  for  evaluating  add-on  solar  systems  for 
use  in  older  buildings.   It  is  conceivable  that  should  the 
results  of  such  programs  demonstrate  a  positive  savings 
through  the  use  of  add-on  equipment,  MIT  would  be  justified 
in  embarking  on  a  limited  program  to  do  the  same.   It  is  not 
envisioned,  however,  that  this  type  of  modular  additivity 
would  ever  serve  as  anything  but  a  supplement  to  the  Central 
Utility  Plant. 

There  is  no  reason  to  exhibit  optimism  about  the  present 
role  of  solar  energy  at  MIT.   Whatever  advantages  there  might 
be  lie  in  the  future.   For  MIT  to  move  now  toward  a  total 
energy  concept  which  includes  solar  measures  is  to  presuppose 
that  generous  incentives  will  be  forthcoming  from  the  federal 
government.   There  is  virtually  no  chance  that  this  will 
happen.   Moreover,  there  are  numerous  questions  relating  to 
urban  sun  rights  I  2 ]  which  have  only  recently  received 
publicity.   At  the  very  least,  decisions  about  the  use  of 
solar  energy  at  MIT  should  wait  until  resolution  of  this 
matter. 

2.2.4   Fuel  Availability 

Inasmuch  as  the  specific  equipment  mix 
which  comprises  a  total  energy  system  presumes  the  use  of  one 
or  perhaps  two  fuels,  it  is  worthwhile  to  examine  briefly  the 
prospects  for  steady  supply  to  the  New  England  area  of  the 
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two  fuels  which  would  most  likely  be  used.   It  is  reasonable 
to  assume  that  those  fuels  in  the  greatest  supply  will 
dominate  the  selection  of  a  power  plant  configuration. 

Because  of  New  England's  unfavorable  location  on  the 
domestic  oil  and  natural  gas  pipelines,  imported  residual  oil 
and  domestic  coal  are  the  principal  fossil  fuels^  which  are 
available  commercially  for  use  in  a  total  energy  scheme  at 
MIT.   Until  1966,  coal  was  the  major  fuel  source  for  electri- 
city generated  in  this  region  of  the  country.   That  year, 
however,  import  controls  were  removed  on  residual  oil.   For 
the  next  seven  years  the  least  expensive  environmentally 
acceptable  fossil  fuel  delivered  to  New  England  was  foreign 
residual  oil;  but,  since  1973,  delivered  oil  prices  have 
exceeded  the  coal  prices  per  unit  of  energy.   When  pollution 
control  costs  for  coal  are  factored  in,  the  two  energy  sources 
appear  equally  attractive  by  most  estimates  [  8  ] .   In  light 
of  the  abundant  resources  of  coal  in  the  U.S.  (90%  of  all 
conventional  energy  reserves) ,  it  is  curious  that  a  larger 
disparity  does  not  exist.   Judging  from  President  Carter's 
expressed  desire  for  industry  to  convert  to  coal,  it  is 
anticipated  that  the  price  differential  will  widen.   Also, 
since  the  world  reserves  of  oil  are  being  depleted  more 
rapidly  than  U.S.  coal  reserves,  foreign  oil  prices  might 
reasonably  be  expected  to  rise  in  the  future  at  a  faster  rate 
than  coal  prices. 

Still,  the  outlook  for  New  England  concerning  supply  of 
coal  is  not  particularly  bright.   The  closest  actively-mined 
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coal  fields  to  the  Boston  area  are  in  southwestern  Pennsyl- 
vania, a  distance  of  approximately  650  miles.   There  are  only 
two  methods  by  which  coal  may  reach  Boston:   rail  and  barge. 
Conceivably,  coal  from  Pennsylvania  could  be  transported  by 
rail  to  a  port  in  New  York,  New  Jersey  or  Connecticut  and 
then  transferred  to  barges  for  the  remainder  of  the  trip. 
More  likely  than  not,  however,  coal  would  be  shipped  by  rail 
[  8 ] .   It  is  interesting  to  note  that  only  one  time  since 
1967  has  a  trainload  of  coal  made  the  trip  north  from 
Pennsylvania.   A  sudden  wide  scale  shift  to  coal  for  the  New 
England  region  would  point  up  at  least  one  pitfall  of  the 
President's  plan  for  coal  conversion:   namely,  that  until 
substantial  improvements  in  the  railroad  track  system  are 
made  in  this  area  of  the  country,  full  scale  delivery  of 
coal  to  potential  industry  users  cannot  be  effected.   The 
present  condition  of  the  rail  system  is  such  that  only  limited 
delivery  schedules  can  be  met  on  a  regular  basis.   What  is  at 
issue  is  the  reclassification  of  the  priority  of  New  England 
railroads  for  rehabilitation  [  8 ]  under  the  provisions  of  the 
Railroad  Reorganization  Act  of  1976. 

The  foregoing  is  not  meant  to  imply  that  coal  is  the 
preferred  choice  of  fuel  for  a  total  energy  system  at  MIT. 
Since  the  Central  Utility  Plant  is  presently  designed  to 
operate  on  fuel  oil,  MIT  could  be  exempt  from  the  requirement 
to  convert  to  coal  if,  in  fact,  cogeneration   is  adopted. 
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2.2.5   Environmental  Impact  of  a  Total  Energy  System 

Should  the  economic  analysis  show  that  the 

on-site  generation  of  electricity  is  the  most  cost  effective 

alternative  for  MIT,  an  environmental  impact  statement  must 

be  prepared  for  the  proposed  plant  configuration.   There  are 

potentially  three  categories  of  air  quality  regulations  to 

which  a  cogeneration  plant  in  Massachusetts  must  conform: 

(1)   federal  and  state  ambient  air  quality  standards 
(AAQS), 


(2)   federal  New  Source  Performance  Standards  (NSPS) , 


and 


(3)  Massachusetts  Air  Pollution  Control  Regulations, 
including  emissions  limitations  and  fuel  quality 
standards  [  9  ] . 

The  federal  ambient  air  quality  standards  for  particu- 
lates, SO  and  NO  ,  were  adopted  as  the  state  standard  by 
Massachusetts.   The  primary  standards  define  the  maximum 
permissible  atmospheric  pollutant  concentrations  which  pro- 
vide for  an  adequate  margin  of  safety  to  the  public.   The 
federal  NSPS  were  promulgated  by  the  Environmental  Protection 
Agency  as  directed  by  the  Clean  Air  Act.   These  standards 
establish  a  maximum  level  of  pollutant  emission  per  unit  of 
heat  input.   NSPS  for  particulates,  S0x  and  NOx  have  been 
promulgated  for  fossil-fuel  fired  steam  generating  units  of 
more  than  250  million  Btu  per  hour  heat  input.   Also  under 
the  provisions  of  the  Clean  Air  Act  is  the  requirement  that 
each  state  adopt  a  plan  which  provides  for  the  implementation, 
maintenance  and  enforcement  of  the  primary  ambient  standards. 


-31- 

Massachusetts  adopted  regulations  on  particulate  emissions 
which  are  actually  more  stringent  than  the  federal  NSPS.   In 
all  cases  of  conflict  between  state  and  federal  regulations, 
the  more  stringent  regulation  is  applicable. 

It  is  quite  obvious  that  the  satisfaction  of  clean  air 
standards  for  the  Boston  area  will  require  considerable 
monitoring  of  pollutant  concentrations.   Allowing  for  possible 
equipment  modifications  to  achieve  acceptable  pollutant 
levels,  the  final  system  cost  will  be  a  function  of  the 
quality  of  fuel  burned. 
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III   STEAM  LOAD  PROFILE  DETERMINATION 

3. 1  Objective 

For  the  purpose  of  total  energy  system  selection, 
a  model  is  required  which  accurately  reflects  the  steam  heat- 
ing and  air  conditioning  loads  at  MIT  over  the  course  of  a 
year.   Ultimately  to  be  incorporated  into  a  computer  program 
which  simulates  the  operation  of  several  total  energy  system 
designs,  the  steam  load  model  must  provide  sufficient 
flexibility  to  allow  a  prediction  of  campus  loads  based  on 
readily  quantifiable  parameters.   Results  will  be  in  the  form 
of  daily  load  profiles  which  describe  the  hourly  variation 
of  campus  steam  demands. 

3.2  Methodology 

From  the  outset  several  factors  were  known  to 
influence  the  Institute's  steam  load.   Chief  among  these 
was  outside  ambient  temperature.   As  steam  space  heating 
is  used  extensively  at  MIT,  the  Central  Utility  Plant  must 
generate  a  steadily  increasing  quantity  of  steam  as  the  out- 
side temperature  drops.   Indeed,  data  records  show  this 
relationship.   Additionally,  wind  velocity,  through  its 
influence  upon  the  heat  transfer  film  coefficient  for  tur- 
bulent flow,  was  known  to  play  a  role  in  increasing  the  heat 
loss  of  buildings.   Not  so  obvious  as  temperature  and  wind 
are  the  effects  of  humidity  and  sun  cover.   Indirect  building 
heat  gains  can  be  attributed  to  both  of  these  factors, 
although  a  precise  determination  of  the  magnitudes  involved 
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is  difficult.   It  was  initially  envisioned  that  the  load 
model  should  account  for  the  effects  of  the  forementioned 
ambient  parameters.   Two  methods  were  considered  by  which 
to  evaluate  the  contribution  of  each  factor. 

The  more  analytic  approach  requires  that  a  detailed 
heat  balance  be  performed  on  each  of  the  campus  buildings. 
By  considering  the  individual  building  construction  and 
determining  a  film  heat  transfer  coefficient  for  windows  and 
exterior  surfaces,  an  overall  heat  transfer  coefficient  may 
be  derived  for  the  walls  and  windows  of  all  MIT  buildings. 
Expressed  in  units  of  Btu/degree-day,  this  information  could 
provide  a  basis  for  evaluating  space  heating  loads  for  any 
particular  degree-day.   Effects  of  heat  loss  through  infiltra- 
tion could  be  estimated  by  the  techniques  outlined  in  the 
ASHRAE  Handbook  of  Fundamentals.   Heat  losses  due  to  wind 
may  be  quantitatively  assessed  by  applying  heat  transfer 
theory  for  forced  convection  over  a  flat  plate  to  the  exterior 
surfaces  of  all  buildings.   Similarly,  heat  gain  through 
window  insolation  can  be  approximated.   As  the  directional 
orientation  of  each  window  is  known,  a  model  could  be  con- 
structed to  yield  building  heat  gain  as  a  function  of  solar 
azimuth  on  a  hourly  basis.   Apart  from  the  space  heating  load, 
the  campus  hot  water  load  could  be  estimated  by  construction 
of  a  hot  water  demand  model,  described  by  a  time  dependent 
usage  function.   In  theory,  therefore,  it  is  possible  to 
analytically  model  the  heating  season  at  MIT,  once  the  ambient 
temperature  and  wind  information  are  available. 
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The  problem  is  somewhat  complicated,  however,  in  the 
warmer  months  of  the  year.   During  the  air  conditioning  season 
steam  load  ceases  to  bear  an  inverse  relationship  to  ambient 
temperature.   Rather,  it  increases  with  ambient  temperature, 
reflecting  the  use  of  the  steam  driven  centrifugal  compressor 
units  of  the  Central  Utility  Plant's  chiller  system.   To 
adequately  deal  with  the  changing  relationship  of  steam  demand 
versus  temperature,  therefore,  an  additional  model  would  be 
needed  to  describe  individual  building  cooling  requirements. 
While  this  does  not,  in  itself,  render  the  load  analysis 
untenable,  the  task  of  steam  load  modeling  using  analytic 
techniques  is  clearly  time  consuming.   Even  more  significant 
is  the  realization  that  it  is  at  best  an  approximation.   In 
trying  to  predict  the  hourly  variation  of  campus  steam  demand, 
there  is  no  guarantee  that  the  magnitudes  so  derived  would 
accurately  reflect  those  which  are,  in  fact,  observed.   An 
alternate  method  is  sought. 

MIT  is  in  the  fortunate  position  of  having  available 
detailed  steam  load  data  for  each  day  of  the  year.   That  is, 
a  recorded  history  exists  of  hourly  steam  demand  as  well  as 
average  wind  and  humidity  conditions  during  the  day.   Using 
this  information  a  variety  of  correlations  may  be  established, 
with  the  result  that  a  load  model  may  be  constructed.   In 
contrast  to  the  analytic  method,  reduction  of  existing  data 
ensures  that  the  magnitudes  of  predicted  steam  loads  are 
representative  of  those  which  would  be  observed  for  any 
particular  degree-day.   Any  disparity  between  what  the  model 
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dictates  and  what  is  actually  perceived  is  thereby  eliminated. 
This  empirical  method  obviates  the  need  to  determine  infil- 
tration heat  losses  or  insolation  heat  gains  since  they  are 
implicity  accounted  for  in  the  historical  steam  data.   Wind 
effects  may  be  evaluated  by  graphing  the  daily  mass  flow  of 
steam  versus  degree-day  for  days  with  and  without  wind.   A 
determination  can  then  be  made  as  to  whether  a  correction 
should  be  applied  to  account  for  an  average  wind  in  the 
Cambridge  area  or  whether  the  wind's  effects  are  negligible 
on  steam  demand.   A  similar  procedure  can  be  employed  to 
determine  the  influence  of  humidity  on  air  conditioning  load 
in  the  summer  months. 

Because  the  purpose  herein  is  not  to  predict  loads  for 
new  buildings  but  rather  to  model  existing  demands,  the 
second  of  the  two  methods  outlined  has  been  chosed  as  the 
more  useful.   Several  assumptions,  however,  are  necessary  to 
permit  modeling  by  this  means. 

3. 3   Assumptions 

Steam  load  is  assumed  to  be  a  function  only  of 
ambient  temperature,  wind  velocity  and  humidity.   The  effect 
of  sun  cover,  which  averages  55%  in  the  Boston  area,  is 
implicitly  accounted  for  in  the  steam  demand.   No  attempt  has 
been  made  to  break  down  the  historical  data  so  days  with 
similar  percentage  sun  cover  are  grouped  together.   It  is 
assumed  that  the  days  chosen  for  data  reduction  comprise  a 
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representative  mix  of  days  with  percentages  of  sun  cover 
typical  of  the  Boston  area. 

It  is  assumed  that  data  from  the  period  January  1976  to 
February  1977  characterizes  present  Institute  steam  loads. 
For  several  years  prior  to  this,  steam  load  decreased  as  a 
result  of  campus  energy  conservation  measures  (see  Section 
3.4).   The  chosen  sampling  period,  however,  represents  a  time 
frame  during  which  loads  have  leveled  out.   Future  steam  load 
growth  will  be  referenced  to  the  above  period  for  use  in  this 
load  model.   A  further  assumption  is  that  the  form  of  the 
daily  load  profiles,  determined  herein,  will  remain  invariant 
with  Institute  growth.   Although  the  magnitude  of  steam 
demand  will  likely  increase  in  future  years,  the  profile 
shapes  will  remain  essentially  the  same.   That  is,  campus 
usage  patterns  will  not  change. 

All  steam  loads  are  to  be  treated  as  one.   Since  the 
major  source  for  building  space  and  hot  water  heat  is  200  psi 
steam,  no  purpose  is  served  by  breaking  down  usage  according 
to  category.   Similarly,  inasmuch  as  the  power  turbines  for 
the  central  air  conditioning  system  use  200  psi  steam  and 
this  represents  simply  another  "load"  on  the  ste<im  system, 
it  need  not  be  segregated. 

3.4   Procedure  For  Data  Collection 

Since  the  oil  embargo  of  1973,  significant  conserva- 
tion programs  have  been  undertaken  at  MIT,  and  they  have 
served  to  lower  the  steam  load  substantially.   During  the 
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period  covering  F/Y  1973  to  F/Y  1976  the  Institute's  steam 
demand,  on  an  annual  basis,  was  reduced  by  26.7%.   At  the 
present  time  there  is  little  room  for  futher  major  reductions 
in  steam  consumption  in  existing  buildings.   In  light  of  this 
information,  it  was  decided  that  for  load  modeling  purposes 
only  the  most  recent  steam  data  should  be  used. 

Information  pertinent  to  steam  loads  at  MIT  is  available 
at  the  Central  Utility  Plant  and  the  offices  of  the  Physical 
Plant.   An  examination  of  load  graphs  at  the  Physical  Plant 
led  to  January  1,  1976  as  the  choice  of  starting  date  for 
data  collection.   The  boiler  operating  logs  at  the  Central 
Utility  Plant  were  used  to  gather  the  daily  total  steam 
generated  for  every  day  during  the  period  January  1,  1976  to 
February  28,  1977.   In  addition,  the  hourly  boiler  steam 
flow  was  noted  for  approximately  30%  of  the  sample  days. 

A  small  percentage  of  the  steam  for  campus  heating  pur- 
poses (notably  for  some  buildings  on  the  East  Campus)  is 
provided  from  Cambridge  Electric  Company.   While  a  daily 
breakdown  of  this  steam  is  not  available,  monthly  totals  are 
available  through  the  accounting  offices  of  the  Physical 
Plant.   A  method  was  required,  however,  to  apportion  the 
monthly  total  of  steam  provided  by  Cambridge  Electric  over 
each  day.   It  was  not  permissible  to  simply  divide  the 
monthly  total  by  the  number  of  days  and  attribute  an  equal 
flow  to  each  day.   This  implies  no  temperature  dependency. 
Accordingly,  the  following  scheme  was  devised: 
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(a)  Over  the  same  monthly  period  for  which  Cambridge 
Electric  steam  data  was  available,  daily  totals  of 
steam  generated  at  the  Central  Utility  Plant  were 
summed. 

(b)  For  each  day  within  this  period  the  decimal  fraction 
of  steam  generated  at  the  Central  Utility  Plant 
relative  to  the  total  from  part  (a)  was  computed. 

(c)  The  decimal  fraction  obtained  in  part  (b)  was 
multiplied  by  the  monthly  total  of  steam  supplied 
by  Cambridge  Electric  to  obtain  the  daily  total  of 
steam  from  Cambridge  Electric. 

Proceeding  on  this  basis  for  the  entire  fourteen  month  sample 
period,  computations  of  daily  steam  furnished  by  Cambridge 
Electric  were  made.   These  were  added  to  the  daily  totals 
from  the  Central  Utility  Plant  generation  to  arrive  at  com- 
posite totals  of  daily  steam  load  for  MIT. 

Although  temperature,  wind  and  humidity  information  is 
available  from  the  operating  logs  at  the  Central  Utility 
Plant,  a  more  reliable  source  was  found  to  be  the  MIT  Depart- 
ment of  Meteorology.   Recorded  using  the  sophisticated  equip- 
ment in  Building  54,  average  wind  and  humidity  data  as  well 
as  average  and  extreme  temperature  information  for  MIT  is 
accessible. 

3. 5   Model  for  Predicting  Daily  Total  Steam  Load 
3.5.1   Preliminary  Data  Reduction 

By  way  of  attempting  to  identify  the 
influence  of  ambient  parameters  (temperature,  wind  and  humi- 
dity) on  daily  steam  load,  a  number  of  graphs  were  constructed 
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to  aid  in  data  interpretation.   Weekdays  were  plotted  separate 
from  weekends/holidays.   The  graphs,  while  not  presented  in 
this  study,  demonstrated  conclusively  that  temperature  and 
wind  significantly  influence  steam  load  while  humidity  has  a 
less  predictable  effect. 

The  first  graph  constructed  was  a  scatter  plot  of  average 
temperature  versus  steam  flow  for  weekdays  with  an  average 
wind  velocity  less  than  10  mph.   Since  the  subject  of  interest 
was  the  determination  of  wind  chill,  the  data  was  limited  to 
days  with  average  temperatures  less  than  65 °F.   An  overlay 
was  made  of  days  with  similar  temperature  and  average  wind 
velocities  greater  than  15  mph.   Although  both  plots  were 
characterized  by  data  groupings  which  suggested  straight  line 
fits,  the  second  graph  was  noticeably  displaced  above  the 
first.   As  the  only  changing  parameter  between  the  two 
groups  of  data  was  the  average  wind  velocity,  it  was  verified 
that  for  days  with  an  average  temperature  less  than  65 °F, 
wind  chill  has  an  augmenting  effect  on  daily  total  steam  flow. 
This  was  to  be  expected. 

According  to  weather  bureau  records,  an  average  wind  of 
approximately  12.5  mph  prevails  in  the  Boston  area.   In  June 
of  1974,  Professor  A.  L.  Hesselschwerdt,  for  a  report  sub- 
mitted to  the  MIT  Physical  Plant,  constructed  a  wind  velocity- 
temperature  correction  chart  [10]  which  is  referenced  to  this 
average  velocity.   It  can  be  used  to  determine  the  equivalent 
temperature  reduction  for  wind  velocities  greater  than  12.5 
mph  as  well  as  equivalent  temperature  elevation  for  velocities 
less  than  the  average. 
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In  order  to  develop  a  means  for  predicting  daily  total 
steam  flow  as  a  function  of  temperature  alone,  the  effects  of 
wind  and  humidity  must  be  allowed  for  as  an  adjustment  to 
ambient  temperature.   Using  Professor  Hesselschwerdt * s 
correction  chart,  this  was  possible  for  wind  velocity.   A 
series  of  scatter  plots  were  initially  made  of  degree-day 
versus  daily  total  steam  flow.   Weekdays  and  weekends/holidays 
were  grouped  separately.   As  a  starting  point,  no  correction 
was  applied  for  wind  velocity.   The  data  showed  considerable 
spread  while  still  suggesting  a  straight  line  data  fit; 
however,  when  the  wind  correction  was  added  for  those  days 
with  average  wind  velocity  other  than  12.5  mph,  the  scatter 
plots  became  significantly  tighter.   To  be  certain,  the  wind 
correction  proved  an  aid  in  data  interpretation  -  to  such  an 
extent,  in  fact,  that  further  refinement  of  the  model  appeared 
possible  using  a  computer. 

Plots  were  also  made  for  days  above  65 °F  of  average 
temperature  versus  steam  flow,  each  point  annotated  with  its 
relative  humidity.   The  intent  here  was  to  determine  what,  if 
any,  correction  should  be  applied  to  average  temperature  to 
account  for  humidity  effects.   It  was  first  thought  that  a 
correction  similar  to  that  developed  for  wind  velocity  might 
result.   The  data,  however,  showed  such  wide  variance  that 
no  correlation  was  possible.   More  specifically,  there  was 
not  even  consistent  evidence  that  higher  humidity  contributes 
to  an  increased  air  conditioning  load.   For  this  reason,  it 
was  decided  to  attempt  a  steam  flow/ temperature  correlation 
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for  days  above  65 °F  with  uncorrected  raw  data.   It  can  be 
argued  that  this  approach  detracts  from  an  otherwise 
rigorous  analysis  of  steam  loads.   The  results,  however, 
indicate  that  it  was  justified  in  that  a  highly  reliable  model 
was  ultimately  developed  which  predicts  steam  demand  over  a 
full  range  of  temperature,  assuming  only  normal  wind  condi- 
tions of  12.5  mph. 

3.5.2   Computer  Analysis 

The  scatter  plots  mentioned  previously 
demonstrated  that  both  weekday  and  weekend  data  would  lend 
themselves  well  to  further  reduction  on  a  computer.   It  was 
initially  thought  that  two  straight  lines  might  best  approxi- 
mate the  data,  one  for  days  with  temperature  less  than  65°F 
and  one  for  days  with  temperature  greater.   The  relatively 
close  grouping  of  data  points,  however,  suggested  that  a 
single  polynomial  curve  fit  might  also  be  possible. 

Using  the  "MIT-SNAP"  program  [11],  daily  steam  flow 
data  refinement  was  accomplished.   Developed  by  the  Sloan 
School  of  Management,  MIT-SNAP  is  an  interactive  data 
analysis  system  for  the  IBM-370  computer.   It  is  designed  to 
perform  basic  statistical  analyses  on  batches  of  data  and 
will  produce  a  least-squares  multiple  regression  of  several 
variables. 

Weekdays  and  weekend/holidays  were  input  as  separate 
groups.   All  ambient  temperatures  were  corrected  to  the  12.5 
mph  wind  velocity  base.   With  temperature  as  the  dependent 
variable  and  daily  total  steam  flow  as  the  independent 
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variable,  a  least-squares  multiple  regression  was  specified 
for  second,  third,  fourth  and  fifth  order  polynomials. 

3.5.3   Presentation  of  Results 

The  best  polynomial  fit  resulted  from  the 
third  order  regression  analysis.   Included  at  the  end  of  this 
chapter  are  the  computer  outputs  for  this  run.   The  equations 
for  the  curve  fits  are: 
for  weekdays: 

S  =  4.1719  x  106  +  32359. 3867T  -  2588.0669T2  +  22.5004T3   (3.1) 

for  weekends/holidays: 

S  =  3.7619  x  106  +  51108. 5156T  -  3081.0232T2  +  26.4316T3   (3.2) 

where 

T  =  average  daily  temperature,  °F 
S  =  daily  total  steam  demand  for  MIT,  lb/day 
For  the  weekday  regression  295  data  points  were  used,  while 
130  were  used  for  weekends  and  holidays. 

A  convenient  feature  of  the  MIT-SNAP  program  is  the  x-y 
plotting  of  all  input  data.   It  can  be  seen  from  Figure  3.10 
and  Figure  3.14  that  a  relatively  tight  grouping  of  points 
exists.   Such  uniformity  is  especially  fortunate  in  view  of 
the  fact  that  wind  velocity  was  the  only  ambient  parameter  for 
which  temperature  was  corrected.   (A  number  in  the  place  of  an 
asterisk  denotes  more  than  one  day  with  that  temperature  and 
steam  flow.)   Despite  the  paucity  of  weekend  data  points  com- 
pared with  weekday,  the  data  trend  is  clear.   The  variance  in 
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daily  total  steam  flow  for  any  specific  temperature  day  can  be 
attributed  to  several  factors.   Foremost  among  these  is  the 
hourly  variation  of  ambient  temperature  during  the  day. 
Although  two  days  may  be  identified  by  the  same  average  tem- 
perature, corrected  for  wind,  one  might  have  excessively  cold 
daytime  temperatures  while  the  other  might  have  warmer  daytime 
temperatures.   Clearly,  the  daily  heating  demand  for  these 
days  could  differ  substantially. 

The  regression  statistics  on  pages  80  and  81  indicate 

2 

the  quality  of  each  curve  fit.   The  high  R  for  both  weekdays 

and  weekends  implies  that  temperature  alone  is  an  outstanding 
predictor  of  daily  total  steam  demand.   In  that  the  F  statis- 
tics are  quite  large,  temperature  is,  indeed,  a  significant 
parameter  in  the  regression  analysis.   The  following  defini- 
tions for  R  and  F  apply: 


R2 


1     —    - 
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k*)2 

R2       1 

1-k 

(3.3) 


F  ■  I=F  TT  (3-4) 


where 


y.  ■  input  data  point 

y  =  fitted  data  point 

y  =  arithmetic  mean  of  all  input  data 

N  =  number  of  data  points 

k  =  degrees  of  freedom 
The  magnitude  of  the  F  statistics  indicates  that  the  variance 
in  steam  demand  explained  by  the  regression  (temperature)  is 
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many  times  greater  than  the  variance  which  is  left  unexplained. 
The  only  mismatch  of  any  consequence  between  the  data  and  what 
is  predicted  by  the  polynomial  approximations  occurs  for  tem- 
peratures less  than  6°F  for  the  weekday  data.   The  error  here 
is  approximately  7%. 

Included  as  Figure  3.12  and  Figure  3.16  are  plots  of  the 
residual  versus  the  fit  for  weekdays  and  weekends  respectively. 
These  are  essentially  displays  of  the  error  between  the  input 
data  and  the  predicted  steam  flow  (fit)  as  a  function  of  the 
magnitude  of  the  predicted  steam  flow.   (The  residual  is 
defined  as  the  difference  between  data  and  fit.)   A  spread  of 
points  dispersed  randomly  about  0.0  on  the  y-axis  indicates 
that  no  other  single  variable  than  temperature  is  necessary  to 
describe  the  variation  of  daily  steam  flow.   This  is  the  case 
for  both  weekdays  and  weekends/holidays  at  MIT. 

A  further  aid  in  data  interpretation  is  provided  by 
Figure  3.13  and  Figure  3.17.   These  show  the  magnitude  of  the 
error  (ordinate)  plotted  against  the  standard  deviation  for  a 
normal  distribution.   On  the  weekday  plot,  for  example,  68.3% 
of  the  data  (one  standard  deviation)  fall  within  an  error 
band  ±  220,000  lb  steam/day.   The  straighter  the  line,  the 
more  even  is  the  error  distribution.   It  is  observed  for  the 
weekend/holiday  data  that  the  error  band  is  not  as  tight  as 
for  weekdays.   This  is  to  be  expected  in  view  of  the  wide 
variation  in  weekend  and  holiday  population  levels  at  MIT 
over  the  course  of  a  year.   Contrasted  to  weekdays  when  MIT 
sees  a  nearly  uniform  number  of  students,  faculty  and 
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management  personnel,  the  weekend  and  holiday  steam  consumption 
is  strongly  a  function  of  the  number  of  students,  alone,  who 
choose  to  remain  on  campus . 

3.6   Seasonal  Variation  in  Daily  Load  Profiles 

From  the  preceding,  a  model  is  now  available  which 
yields  daily  total  steam  flow  as  a  function  of  outside  average 
temperature.   It  may  be  applied  for  any  day  of  the  year. 
Numbers  which  reflect  daily  total  steam  consumption  levels, 
however,  are  themselves  of  little  practical  use.   It  remains 
to  develop  a  method  whereby  this  24  hour  total  may  be  appor- 
tioned over  each  hour  of  the  day.   Indeed,  it  is  this  hourly 
fluctuation  in  steam  demand  which  will  ultimately  govern  the 
selection  of  equipment  for  a  specific  total  energy  system 
alternative. 

3.6.1   Procedure 

Blocks  of  days  were  chosen  from  the  fourteen 
month  sampling  period  as  being  representative  of  winter,  fall, 
summer  and  spring.   Three  weeks  of  data  were  obtained  for  each 
season.   For  each  day  the  hourly  steam  demand  as  transcribed 
from  the  Central  Utility  Plant  operating  log  was  normalized 
with  respect  to  the  hourly  average  for  that  day.   Hereafter 
referred  to  as  the  "hourly  load  factor"  method,  this  simplifi- 
cation provided  an  efficient  means  of  data  quantification. 

The  choice  of  grouping  days  by  season  was  made  somewhat 
arbitrarily.  It  was  hypothesized  that  daily  steam  load  pro- 
files might  follow  a  seasonal  pattern  and  that  attempts  at 
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load  modeling  should  initially  concentrate  on  defining 
repeatable  profiles.   Proceeding  on  this  basis,  crude  graphs 
of  hourly  load  factor  versus  hour  of  day  were  made  for  several 
days  within  each  season.   These  served  to  verify  that,  in  fact, 
seasonal  profiles  did  exist  for  both  weekdays  and  weekends. 
With  the  foundation  thus  set  for  further  refinement  of  the 
data,  more  sophisticated  techniques  were  employed  to  determine 
the  representative  profiles  for  each  season. 

3.6.2   Attempts  at  Polynomial  Approximations 

For  the  purpose  of  modeling  it  was  envisioned 
that  an  analytic  expression  which  described  the  variation  in 
hourly  load  factor  might  prove  instrumental.   To  this  end, 
the  MIT-SNAP  program  was  utilized  in  a  fashion  similar  to  that 
described  in  Section  3.5.3. 

For  the  first  computer  run,  several  weeks  of  winter  week- 
day data  were  used  in  the  regression  program.   Second,  third, 
fourth  and  fifth  order  multiple  regressions  were  specified. 
While  MIT-SNAP  did  provide  a  graphical  plot  of  all  input  data, 
the  analytical  results  were  less  than  satisfactory.   Because 
the  program  incorporates  a  matrix  inversion  feature,  a  high 
degree  of  colinearity  was  found  to  exist  among  the  coefficients 
of  polynomial  approximations  of  order  three  or  higher.   Con- 
sequently, the  equations  generated  by  MIT-SNAP  were  not 
equations  which  described  the  input  data.   A  second  program 
was,  therefore,  developed  which  relied  upon  an  IMSL  library 
subroutine  (LSFIT)  to  perform  a  least  squares  regression. 
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More  positive  results  were  obtained  from  the  use  of  the 
LSFIT  subroutine.   With  hour  of  the  day  specified  as  the 
independent  variable,  the  regression  was  carried  out  for  each 
seasonal  group  of  data.   The  higher  order  polynomial  equations 
(fourth  and  fifth  order)  for  the  approximating  curves  matched 
the  data  well  in  some  cases.   There  was  a  consistent  disparity, 
however,  between  the  fitted  curve  and  data  for  the  hours  8:00 
P.M.  through  12:00  P.M.   Also,  in  several  instances,  the  peak 
magnitude,  as  predicted  by  the  polynomial,  was  substantially 
less  than  the  data  would  seem  to  indicate  that  it  should  be. 
Although  repeated  attempts  were  made  to  manipulate  the 
approximating  equations  so  as  to  obtain  more  exact  fits,  it 
became  all  too  obvious  that  a  purely  analytic  means  of  modeling 
would  not  be  possible. 

3.6.3   Hourly  Scale  Factor  Adjustment 

With  the  exception  of  only  several  hours  in 
any  one  seasonal  profile,  the  hourly  scale  factors  derived 
using  the  polynomial  approximation  techniques  were  most 
representative  of  the  input  data.   It  was  found  that  by 
adjusting  the  magnitude  on  some  of  the  scale  factors  so  as  to 
more  accurately  reflect  the  seasonal  trend,  the  remaining 
disparity  between  input  data  and  the  curve  fit  could  be 
eliminated.   This  procedure  was  attempted  for  the  purpose  of 
determining  what  increase  in  correlation  coefficient  could  be 

achieved  over  that  resulting  from  the  polynomial  approximation 

2 
alone.   The  higher  the  R  ,  the  more  closely  a  data  group  is 

described  by  a  specific  "curve  fit" . 
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For  the  weekday  and  weekend/holiday  steam  profiles, 

hourly  scale  factor  adjustments  were  made  and  new  values  of 

2 
R  were  computed.   An  iterative  process,  scale  factor 

corrections  were  made  so  as  to  achieve  the  highest  possible 

R  ,  i.e.,  minimum  error  sum  of  squares,  consistent  with  the 

2 

general  trend  of  the  data  input.   Increases  of  R  in  the  range 

of  .05  to  .20  resulted.   Although  the  profiles,  after  the 
scale  factor  adjustment,  were  no  longer  smooth  as  the  poly- 
nomial approximations  would  have  them,  the  resulting  fit  was 
in  each  case  highly  consistent  with  the  data. 


3.6.4   Weekday  Results 

Included  as  Figures  3.1  through  3.5  are  the 
final  steam  load  daily  profiles  for  MIT.   Each  represents  a 
composite  profile  inasmuch  as  it  reflects  a  balance  between 
polynomial  approximation  techniques  and  optimization  efforts 
to  ensure  a  minimum  error  sum  of  squares  between  data  and  fit. 
For  any  particular  temperature  day  the  hourly  distribution  of 
steam  demand  may  be  predicted  as  follows: 

(a)  Compute  the  daily  total  steam  demand  as  a  function 
of  temperature  from  equations  3.1  and  3.2  in  Section 
3.5.3. 

(b)  Divide  the  daily  total  from  part  (a)  by  24  to 
arrive  at  the  average  hourly  demand. 

(c)  Multiply  the  average  hourly  demand  by  the  respective 
hourly  load  factor  to  determine  the  "predicted" 
steam  load  for  that  hour. 
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As  the  seasonal  profile  differences  are  distinct,  several 
observations  can  be  made  regarding  the  weekday  graphs.   Of 
interest  is  the  variation  in  magnitude  of  peak  hourly  load 
factor  over  the  course  of  a  year.   Although  the  winter  peak 
is  less  than  that  for  any  other  season,  the  actual  magnitude 
of  hourly  demand  is  the  highest.   Thus,  the  profiles,  by 
themselves,  provide  no  absolute  information  on  loads;  rather 
they  describe  only  the  relative  fluctuation  in  demand  with 
respect  to  a  daily  average.   Note  also  the  steep  daily  peak 
around  10:00  A.M.  for  both  summer  and  fall  weekdays.   Largely 
due  to  the  campus  air  conditioning  load,  these  peaks  reflect 
the  initial  daily  surge  in  cooling  demand  which  accompanies 
the  arrival  of  the  MIT  community  in  the  morning  hours.   As 
outside  doors  are  opened  and  buildings  which  were  closed 
during  the  evening  assume  their  normal  occupancy  levels,  the 
Central  Utility  Plant's  chiller  system  requires  an  increasing 
amount  of  steam  to  operate  its  turbine  driven  compressors. 
Once  the  load  stabilizes,  the  system  responds  by  a  somewhat 
reduced  steam  demand  for  the  remainder  of  the  day. 

The  fact  that  only  one  profile  each  exists  for  summer  and 
fall  weekdays  does  not  imply  that  this  profile  alone  is  pre- 
cisely repeated  day  after  day.   Certainly,  there  are  variations 
in  all  of  the  seasons.   Although  the  time  of  occurrence  of 
daily  peaks  is  strikingly  similar  within  any  one  season,  the 
peak  magnitudes  vary.   The  goal  here,  however,  has  been  that 
of  identifying  the  predominant  profile (s). 
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Both  winter  and  spring  weekdays  were  characterized  by 
two  distinct  demand  patterns.   The  normal  pattern  results 
when  steam  demand  is  roughly  proportional  to  building  usage. 
That  is,  beginning  at  approximately  9:00  A.M.  and  continuing 
until  6:00  P.M.,  steam  load  is  greater  than  any  other  time  of 
the  day.   A  review  of  hourly  ambient  temperature  fluctuations 
were  made  for  days  in  this  category.   It  showed  that  tempera- 
ture remains  relatively  constant  during  the  daylight  portions 
of  such  days.   The  extreme  pattern  is  characterized  by 
excessive  hourly  load  factors  in  the  early  morning  hours. 
After  approximately  10:00  A.M.  the  scale  factor  variation 
is  similar  to  that  for  the  "normal"  day.   An  explanation  for 
this  anomoly  is  that  during  cold  weather  periods  it  is  not 
uncommon  for  certain  days  to  exhibit  unusual  temperature 
variations.   More  precisely,  when  the  coldest  temperature 
occurs  in  the  early  morning  hours  and  the  temperature  then 
increases  during  the  remainder  of  the  day,  it  is  likely  that 
such  an  extreme  load  profile  will  result.   Approximately  30% 
of  the  winter  and  spring  days  sampled  during  January  1976 
to  February  197  7  exhibited  the  extreme  profile. 

3.6.5   Weekend/Holiday  Results 

Figures  3.6  through  3.9  depict  the  hourly 
load  fluctuation  of  steam  for  weekends  in  each  of  the  four 
seasons.   Although  not  presented  in  this  study,  the  data  from 
which  these  graphs  were  constructed  showed  less  sharply  defined 
profiles  than  for  the  weekdays.   This  is  partially  attributable 
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to  the  fewer  number  of  sample  days  available  for  data  reduc- 
tion.  Additionally,  weekends  and  holidays  have  demand 
patterns  which  are  inherently  less  likely  to  be  repeated  week 
after  week. 

3.7  Steam  Load  Profile  Summary 

Presented  in  Tables  3.1  and  3.2  are  listings  of 
the  hourly  load  factors  for  each  daily  steam  profile.   As  an 
aid  in  determining  the  magnitude  of  hourly  Institute  steam 
demand  for  any  particular  season,  they  are  included  as  a 
supplement  to  the  graphical  representations.   Table  3.3  pro- 
vides a  numerical  listing  of  the  daily  steam  consumption/ 
temperature  information  which  is  reflected  by  Figure  3.10. 
The  same  information  for  weekends/holidays  (Figure  3.14)  is 
shown  in  Table  3.4. 
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Hour 

of 

Day 

Normal 
Winter 
Weekday 

Extreme 
Winter  & 
Spring 
Weekday 

Normal 
Spring 
Weekday 

Normal 
Summer 
Weekday 

Normal 

Fall 

Weekday 

1 

.915 

1.057 

.910 

.840 

.820 

2 

.920 

1.050 

.890 

.825 

.780 

3 

.935 

1.047 

.880 

.810 

.760 

4 

.925 

1.035 

.880 

.790 

.750 

5 

.915 

1.026 

.900 

.780 

.740 

6 

.930 

1.031 

.910 

.795 

.750 

7 

.970 

1.067 

.980 

.840 

.810 

8 

1.030 

1.164 

1.150 

.900 

.910 

9 

1.120 

1.185 

1.210 

1.060 

1.060 

10 

1.150 

1.110 

1.160 

1.280 

1.150 

11 

1.130 

I.060 

1.100 

1.190 

1.270 

12 

1.090 

1.000 

1.080 

1.210 

1.240 

13 

1.080 

.998 

1.060 

1.200 

1.290 

14 

1.050 

.966 

1.040 

1.210 

1.250 

15 

1.035 

.954 

1.000 

1.180 

1.230 

16 

1.020 

.951 

1.000 

1.140 

1.190 

17 

1.010 

.9^0 

1.030 

1.170 

1.170 

18 

1.000 

.937 

1.020 

1.110 

1.150 

19 

.990 

.922 

.980 

1.060 

1.110 

20 

.990 

.938 

.970 

1.000 

1.020 

21 

.970 

.928 

.970 

.960 

.950 

22 

.960 

.897 

.970 

.920 

.900 

23 

.940 

.882 

.960 

.875 

.860 

24 

.925 

.855 

.950 

.855 

.840 

Table  3.1  -  Weekday  Steam  Profile  Hourly  Load  Factors 
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Hour 
of 

Day 

Normal 
Winter 
Weekend  & 
Holiday 

Normal 
Spring 
Weekend  & 
Holiday 

Normal 
Summer 
Weekend  & 
Holiday 

Normal 
Fall 

Weekend  & 
Holiday 

1 

.967 

.950 

.920 

.940 

2 

.964 

.970 

.905 

•  932 

3 

.960 

1.010 

.892 

.924 

k 

.957 

1.010 

.888 

.920 

5 

.955 

1.020 

.888 

.916 

6 

.953 

1.020 

.900 

.918 

7 

1.010 

1.030 

.914 

.920 

8 

1.050 

1.050 

.930 

.935 

9 

1.060 

1.070 

.975 

.980 

10 

1.075 

1.090 

1.011 

1.060 

11 

1.065 

1.060 

1.035 

1.120 

12 

1.020 

1.030 

1.057 

1.150 

13 

1.030 

1.010 

1.109 

1.130 

14 

1.030 

.960 

1.144 

1.090 

15 

1.030 

.920 

1.134 

1.075 

16 

1.015 

.900 

1.154 

1.070 

17 

1.000 

.990 

1.140 

1.090 

18 

1.000 

1.040 

1.123 

1.050 

19 

.980 

1.020 

1.151 

.990 

20 

.980 

.980 

1.060 

•  975 

21 

.980 

.980 

1.000 

.960 

22 

.978 

.970 

.945 

.955 

23 

.972 

.960 

.932 

.950 

24 

.969 

.960 

.900 

.950 

Table  3.2  -  Weekend  Steam  Profile  Hourly  Load  Factors 


■63- 


t~  —  oooooo 

<  Id  U  U  U  Ul  U 

O  X  i/%  o  <*  a 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


oooooo 


HO       o 


l/>  O 

-J  * 


oooooooooooo 


t-      p-  c  tn  o  «/>  ip  «- 
«  c 


j\  r-  i/>  »  m  o  o 


ooooocioocoooooooooooooooooooooooo 


«o  <r   J  C  3  <*•  c. 


MrlMa    l»l  O  "1  r*  «-   » 


Klf'rMfinrMl 


iOifloc~o<»ir-»w»cc»ir<'>o«j>ciNir>or-of>io«N<^ 


a> 

H 

EH 


(OMMflO 


■64- 


oooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

u  (•.  i*i  ui  iii  i<:  ii.  u  (*)  u  w  u  u  i.i  i 


u;  l'  f«:  u  i*:  w  M  l.i 


©OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO© 


Oi/»p*o«^>o*n'^^i/i^>ooo^ 


t*    &  ±3    a    »«<£0^t^' 


^o*-^^(Nf>*tra>*n 


»-  o  r*  o  o  r*  or  o  o  . 


^ir^a^ 


4tAtT   4    4 


o  *>  r-  r*  u"»  ^  -rf>  ■*  j^»rw-4ir<iP>£>vo^ir>rtf>-a  i 


-65- 


OOOOOOOOOOOOi 


oooooooooooooooooooooooooooooooooooooooo 


« .   *:   «    J   J"  r-    J    -    cr»  <T»  r-   o    .*,  ^»  ' 


o   n   i  ^    m  o>  n  i'.  Mr.  o  ff  ff  3   ^r^o>*r*n<r»naoo 
*»JHi-5r<-j,.1,/lO'<MN«e3*l(lON03 


o  o  o  o  o  o 


oc  oooooooooooooooo 


c O  O     i  o 

»  *  <^  o  « 


o  o  o  o  r,  o 


ooooooooooooooooooooooooooooo 


ooooooooo 


oooooocooooooo 


-66- 


OOOQOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

id  id  u:  u.  id  ui  td  u  uj  w  w  m  MM  id  u»  u  id  u  id  id  id  co  ui  w  ui  ik)  r.>  uj  uj  id  ui  w  1*1  us  u  u>  ui  ui  td  ui  K  td  ui  td  uUbiUUWidMUUU 

^•-^^^^^•-^^t-r-«-r-^^^«-^^»-*-»-»-^^»-^^»-r-^^^f\r^f>ir^r^r>iP>iP4rsicsr^f>iiNf>iMWCNr*iNf»i#n*n 
oooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

<*  e-  *>  en  i%  o»  ^  »»»•*«  m  r-  «  *j--»nr*a><oo<»f«o«-»©<**>«-»-e^*>'^<»oor-o^  .*  <n  ©»  o«  *  »  «N*h*o^coh*inMj 


■67- 


oooooooooooooooooooooooooooooooooooooooooooooooooooooo  O'O 
ooooooo  ocoooc  ooooooooooooooooooooooooooooooooooooooooooo 


■68- 


ooooeooooooooooooooo 

UUUWUIdU   UUUUUIilUUUUUUU 

ootftfcntfiho^tf^r*  «  ©  o*  o  r*  »-  *^ 
oooooooooooeoooooooo 


r> 


«»•  Old 


-69- 


»-  fi^ooooeoooooooooooooooooooooooooooooooooooooooooooooo 

a  <UUUI-IUUUUMI.lU.UI»b!UUIUUbllilUUUIItll>lUUUJI>IUUICUU]UU[Ob]UUUIUI>IUUUlill.lluUU 


ooooooooooooooooooooooooooooooooooooooooooooooooooo 


o  o 

tfi  m 

MO 


B.   U         C 
O  O        CI 


-3- 
cd 

Eh 


IflOlf  INi 


oMNMh  oo«r-*^of-»/!«*»of>*o*4r»o< 


ooinr-o<»>>ooo<oaooo<j>o»0 


70 


oooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


UWUUU.UUUU 

»c  «  *.  «  o  •-  «.  o  n 


UWldbltilblUlUUIUUUtiJUUUUUllOblUUUUlMUUbll 

^  u"   juior»^(Nr*oo^io»3>ocDooiii,*f>it"'in«*ir»n»»-mooi 
»r»aa<in«>ON»-3iw«rMo>»»»oeDmif)»-<No»-oo«)r»' 

r«ifN<N<Me«3«»!r*r-r».Tir».r»c>^-ooooifM«-or>rwui»>r'oi 


t-=»fna>«-i-»-oie>or-»mr»<oo 


OOOOOOOOOOi 


ooooooocoooooooooooooooooooooooooooooooo 


H 


o  o  o  o  o  i 
•  «  <r>  in  in  < 


©  O  O  O  O  O  I 

m  r-  <*  «"»  »  O. 


leeeeeeoa 


>/"•  o  o  Jl  vC  or    on  i 
i  <o  a  r*  r-  i 


iiu>»nr  3  irotNicor-oo< 


-71- 


O  ©  O  O  O  I 

III  M  1.1  uj  ul  I 


OOOfJOOOOOOOOOOO 

I  1.1  U  I.]  HI  liM  Ul  1.1  U  U  U  Ul  tl  U  III 


^^c^r^aj^-i/wrui. 


»h  nifino 


OOOOOOOOOOOOOOOOOOOOOOO 


/i  Ji    I   r,   4  .-   J   ai  lA  f  i  ^  c    j'   m  u  Jl  1  <  O  O  I 


-72- 


•  o»       » 

<M 

•  •    •  (M 

04  »  *M 

rtoit  { 

CM  Ol  •  < 

rt  run 

<M  »  N 

•  in  in 


a 
u 

0) 

p. 
e 

a> 

EH 

-p 

C 
0) 


I 


M*         CM  • 


«  • 


s 
s 

d> 
•p 
en 

& 

4) 

Eh 
M 
S 

I 

O 


(XBp/qx)   uondamsuoo  treats  ^TPSCI   TB^oji 


73- 


u 

o  CO 
0) 

0>  P« 

N  I 

•P  En 

CO 

fc  CO 

0)  3 

Ph  to 

S  h 

0)  0) 

Eh  > 


•P 
W 

•3 

•P 

5? 

0) 


..  * 

•r  to  9 

S  K  | 

?  T  -H 


(Xep/qx)   uo-j^dumsuoo  oreads  pa^oxpaad 


-74- 


•    •  • 


•         •  Mr. 


•         •    • 


♦  m 


•  «  • 


•  •  •  r«  *   • 


•  •        •       «■» 

«    •        »  •  • 

•  •     •    <N  •  <"> 


s? 


W 
>> 

i 

< 

c 
o 

•H 
09 
0) 
0) 

fe 

s? 

•a 

M 
0) 
0) 

•S 

U 

o 

<H 

•H 

ft, 

0} 
(0 

> 
m 

3 

•d 

•h 
to 
a> 

oc 

<H 

o 
o 

ft. 
I 
CM 


a> 

(4 

•H 

ft* 


puBuiaa  oreads  XBpx©«M  Joj  sinnpisay 


-75- 


o 

CO 
•H 

•H 

(0 

+» 

>» 

at 

H 

•H 

flj 

> 

c 

a> 

^ 

Q 

>> 

TJ 

03 

is 

T3 

0) 

§ 

a> 

3 

+» 

m 

U 

o 

«M 

c 

o 

.H 

p 

3 

.O 

•H 

U 

•P 

OQ 

•H 

Q 

U 

© 
e 
o 

2 

o 

a 

<D 

> 

•H 

„ 

■P 

O          M 

O          ■ 

o      c 

1 

©     •» 

o 

©       «c 

z 

IT         VI 

u 

o 

H 

p 

M         V) 

o 

B       •■» 

H 

M         £ 

(U 

U        D 

1 

K 

<*> 

« 

«H 

M 

• 

r\ 

<D 

V»         VI 

M         M 

I 

i        i 

(top/qi)   uoT^doinsuoo  urea^s 
P®**Td  PU*  ^i^Q  ueaMq.ag  aojjg 


-76- 


o  0) 


*     «     •     <N 

«     • 


K    P 


(tep/qX)   uo-f^domsuoo  urea^s  ^IT1^   T**oj 


-77- 


a> 
U 

^       -P 
pin       rt 

>-^         0) 

ft 

<d     e 

u  <D 

3  E-< 

•P 

as  oo 

U  3 

a>  to 

ft  fc 

S  a> 

a  > 

EH 


•  « 


..  a> 

£  *» 

«  +> 

°  .h 

5  * 

u  I 
n 


o 


(£«p/qi)    uoi^dumsuoo  aree^s  pe^oxpsjd 


■78- 


»         • 


•         • 


•         « 


o 

o 
o 

o 

& 

(0 
>> 
H 

If) 

T* 

at 

> 

5 

H 

C 

o 

T3 

•H 

o 

* 

§ 

(0 

to 

« 

6 

a> 

© 
© 

tt) 

Q 

& 

o 

0) 

a 

6 

3 

tt 

0) 

•a 

+» 

c 

to 

•o 

a> 

VO 

c 

0) 

o 

CD 

3: 

to 

.* 

o 

0) 

u 

© 

0) 

o 

© 

3 

<N 

« 

XJ 

■P 

(D 

•H 

■P 

fe 

+> 

•H 

CO 

Pn 

3 

CO 

u 

M 

© 
© 
© 
o 

0) 

> 

CO 

3 

o      o 
o      e 
©      © 

©      © 

CO 

©      © 

a> 

>© 

©      o 
©      © 

oc. 

o 

©      © 

J 

in      id 

<H 

o 

o 

o 

© 
o 

H        £ 

■p 

o 

K          ■ 

D         O 

H 

hi         M 

a« 

-»         .J 

«         •* 

i 

in       v> 

..    ..    u 

i/i      m      3 


pmauiea  urea^s  pua3(adM  aoj  siBnpfseH 


-79- 


■n 

4f»  • 


o 

C         (0 

O           -H 

© 

•H           W 

c 

P           >> 

£ 

«J         rH 

>     1 

0)         < 

Q 

T3 

•d       C 

fc       a> 

-a       a) 

C         0) 

© 

rt        5 

2 

•P 

• 

W         fc 

"* 

o 

<n 

c 

o 

•H 

■p 

3 

.o 

•H 

© 

b 

o 

-P 

. 

03 

© 

•H 

Q 

h 

§    ? 

o      fc 

■:   ^ 

©  w 

d> 

o 

> 

o 

•H 

© 

..     -P 

CM 

©     R     S 

©SB 

°      B      H 

O         ft)         O 

o     «i     2 

©     u 

«      ■»     <H 

o 

o 

m      in      a 

© 

B         •■»       iH 

c 

«  c  * 

a 

->      m 

I 

«      •*       1 

t)        B          ' 

10         © 

.     «^ 

S    *"• 

m      m     F^i 

in      b 

n      •« 

sua) 

S  3) 


(^Bp/qx)    uo-jq.duinsuoo  urea^s 


-80- 


«•»       sr  sc  if 

K        O  «N  f- 

O  VC  d 

if.  in  • 


■p 
co 

& 

0) 

a> 

(4 
O 

<H 

CO 

o 


r-        cr-       o>       t- 
©       o       ©       t- 


>       in  ►>.        o       «*■        IS 

M      <e  «       «n      - 

CO  ">       ©        » 


6       in 

»->       » 
to       © 


m       «t       «t       «- 


o  HO 

♦  M        ♦ 

*-  MM 

»  O        O- 


J        an       —        - 


o       o 


©      «r      o> 


M         •* 

m       *a 

2    S 


m-      m.  m  w 

M         *  ft)  M  ft] 

6o  V  C.  -J 

U  A. 


9  •» 

M 

ft. 

O 


-81 


•-       m      r» 


r-       «c       m 


6 

id 

0) 

■p 

M 

w 

M 

as 

•a 

s 

o 

M 

0) 

H 

0) 

s 

s 

m 

u 

O 

o 

* 

O 

(0 

•a 

M 

o 

•H 
•P 

• 

(0 

M 

•H 

M 

\o 

cd 

(9 

r» 

0 

CD 

+» 

10 

r» 

O^ 

cr> 

^ 

w 

►J 

o 

o 

© 

i 

N 

in 

c 

« 

o 

OS 

•H 

H 

CO 

cn      in 

* 

r* 

a 

CO 

o> 

>«      tn 

<N 

o 

e»> 

H 

«r> 

a> 

IT. 

c- 

<D 

■ 

or 

a> 

n 

o 

O 

K 

r* 

o 
i 

o 

<N 

■ 

5b 

m 

u> 

er> 

a) 

H 

•" 

i 

ce 

•-      r-      «o      m 


O 


•-•      p      ©      *" 

j.    8    o 


n  o 

»■  n 


•J  «■ 


a> 


-82- 
IV    ELECTRICAL  LOAD  PROFILE  DETERMINATION 

4.1  Objective 

The  motivation  behind  a  determination  of  specific 
electrical  usage  profiles  at  MIT  is  the  accurate  modeling  of 
representative  electrical  loads  for  use  in  total  energy  system 
selection.   As  with  the  steam  demand  model,  the  ultimate  goal 
of  this  effort  is  the  computer  simulation  of  actual  operating 
loads  at  MIT.   Daily  load  profiles  are  required  which  describe 
the  hourly  variation  of  campus  electrical  demand. 

4.2  Methodology 

Unlike  steam  demand,  electricity  consumption  at  MIT 
is  not  largely  a  function  of  ambient  parameters.   Operation  of 
lights  and  major  venilation  equipment  is  virtually  independent 
of  temperature,  sun  cover  and  wind  speed.   Although  it  is  true 
that  during  warmer  weather  a  decided  load  increase  can  be 
observed,  reflecting  the  wide  spread  use  of  fans  and  space  air 
conditioning  units,  an  accurate  prediction  of  the  magnitudes 
involved  is  most  difficult. 

What  ultimately  governs  electricity  consumption  is  the 
usage  factor  of  each  individual  campus  building.   The  deter- 
mination, therefore,  of  a  method  for  predicting  hourly  load 
variation  is  dependent,  at  the  very  least,  upon  the  successful 
modeling  of  building  occupancy  levels.   A  truly  accurate  load 
assessment  would  require  the  detailing  of  specific  usage 
patterns  for  any  given  level  of  occupancy.   Needless  to  say, 
such  an  undertaking  could  well  prove  exhausting  and  possibly 
fruitless. 
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An  alternative  means  of  load  estimation  centers  on  the 
reduction  of  already  existing  electrical  data  for  the  purpose 
of  determining  what,  if  any,  correlations  may  be  established. 
To  this  end,  information  from  the  electrical  logs  at  the  MIT 
Physical  Plant  was  used  for  the  period  January  1,  1976  to 
February  28,  1977.   In  that  the  vast  majority  of  campus 
energy  conservation  measures  were  implemented  prior  to  this 
time  frame  (reduced  lighting  levels,  equipment  cycling,  etc.), 
the  data  represents  a  stable  base  period. 

4. 3   Model  for  Predicting  Daily  Total  Kilowatt  Demand 

Figure  4.1  is  a  scatter  plot  of  daily  total  kilowatt 
load  versus  ambient  temperature  for  weekdays .   Figure  4.2  is 
the  weekend  plot.   It  can  be  seen  that  a  considerable  varia- 
tion in  electricity  consumption  exists  for  any  particular  band 
of  temperature.   One  reason  for  this  spread  lies  in  the  fact 
that  usage  patterns  at  MIT  are  not  the  same  from  one  day  to  the 
next.   It  is  true  that  a  fundamental  electrical  load  prevails 
each  weekday  and  weekend.   Corridor  lighting  and  dormitory 
ventilation,  for  example,  are  forms  of  demand  which  remain 
relatively  constant  each  day.   On  top  of  this  category  of 
permanent  load,  however,  are  superimposed  a  variety  of 
transient  loads.   For  instance,  classes  which  meet  Monday, 
Wednesday  and  Friday  usually  do  not  meet  Tuesday  and  Thursday. 
Evening  lectures  and  seminars  are  scheduled  unpredictably  and 
at  unevenly  spaced  intervals.   It  is  clear,  therefore,  that 
whatever  model  is  developed  must  address  the  problem  of 
spurious  usage  patterns. 
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The  data  indicates  that  while  daily  total  kilowatt  load 
cannot  be  reliably  predicted  simply  on  the  basis  of  tempera- 
ture, it  is  bracketed  by  certain  upper  and  lower  bounds  of 
consumption.   Days  with  temperatures  less  than  60 °F  show 
essentially  the  same  upper  and  lower  bounds  while  the  kilowatt 
load  for  warmer  days  tends  upward.   Although  there  is  no  one 
explanation  of  why  daily  consumption  varies  so  widely,  the 
fact  that  it  does  vary  between  definable  extremes  suggests 
that  a  scheme  might  be  devised  to  predict  daily  total  kilowatt 
load  based  on  historical  distribution  patterns.   It  is  con- 
ceivable that  a  random  selection  procedure  might  be  employed 
for  the  purpose  of  daily  total  kilowatt  demand  modeling. 

One  alternative  would  be,  for  days  with  temperature  less 
than  60 °F,  the  random  assignment  of  a  specific  daily  kilowatt 
total  so  as  to  fall  within  the  well-defined  upper  and  lower 
bounds  of  Figures  4.1  and  4.2.   For  temperatures  greater  than 
60 °F,  assignment  would  center  around  a  monotonically  increas- 
ing function  of  temperature,  also  incorporating  some  specified 
bandwidth.   While  certainly  a  simplified  approach,  this  method 
implicitly  accounts  for  the  "unknown"  factors  which  cause  the 
day  to  day  variation  in  kilowatt  load.   Its  disadvantage  lies 
in  the  fact  that  nothing  can  be  said  regarding  how  representa- 
tive the  projected  loads  are  of  actual  consumption  data.   A 
means  of  ensuring  a  more  representative  spread  in  daily  kilo- 
watt totals  would  be  to  use  a  well  established  statistical 
distribution  in  making  load  assignments.   By  imposing  the 
tenets  of  the  Central  Gaussian  Theorem  it  can  be  inferred  that 
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with  increasing  numbers  of  data  points,  the  kilowatt  bandwidths 
depicted  in  Figures  4.1  and  4.2  would  mirror  a  normal  distri- 
bution about  some  mean  daily  kilowatt  load.   If  it  were 
required,  therefore,  to  estimate  the  daily  total  kilowatt 
demand  for  a  finite  number  of  days  with  temperatures  less  than 
60 °F,  load  assignments  would  be  made  purely  on  the  basis  of 
proximity  to  the  "mean".   68.3%  of  the  days  would  have  daily 
total  kilowatt  loads  falling  within  one  standard  deviation  of 
the  mean  as  determined  by  the  Gaussian  distribution  within 
the  respective  bandwidth.   Still  a  third  method  would  be  to 
assign  kilowatt  loads  to  days  using  the  same  proportionate 
distribution  as  that  which  characterizes  the  sampling  data. 

Because  the  load  model  is  to  be  used  specifically  for 
simulating  "typical"  consumption  patterns,  it  was  decided  that 
the  third  method  above  should  be  used.   It  is  the  most  con- 
servative of  the  three  approaches  in  that  it  presumes  only 
that  the  specific  mix  of  daily  kilowatt  totals  for  January 
1976  to  February  1977  is  representative  of  the  variation  for 
any  particular  time  span. 

Figures  4.3  through  4.11  show  the  distribution  of  daily 
total  kilowatt  load  for  the  sampling  data.   For  temperatures 
greater  than  60 °F  the  weekday  demand  tends  upward  with  a  slope 
of  18  35  kilowatts/degree,  while  the  weekends  show  a  slope  of 
2600  kilowatts/degree  above  70 °F.   These  percentage  distri- 
butions of  kilowatt  load  may  be  used  to  predict  a  spread  of 
representative  electrical  consumption  totals  for  any  number 
of  days,  providing  only  that  the  temperature  band  is  known. 
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Figure  ^.10 


Kilowatt  Hours 


Distribution  of  Daily  Total  Kilowatt 
Demand  at  MIT  for  Weekends  with 
70°F<  Temperature  £  75°F. 
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4.4  Daily  Load  Profile  Determination 

Data  for  every  second  day  during  the  fourteen  month 
sample  period  was  used  in  the  analysis  of  typical  daily 
electrical  demand  profiles.   Following  the  same  approach  as 
for  the  steam  load  model,  an  average  hourly  kilowatt  demand 
was  determined  for  each  day  and  used  to  compute  hourly  kilo- 
watt load  factors. 

Contrasted  to  the  steam  profiles,  very  little  seasonal 
similarity  was  found  to  exist  for  the  daily  electrical  load 
profiles.   Attempts  were  made  to  identify  repeatable  character- 
istics of  the  data  in  the  hope  of  establishing  a  distinctive 
trend.   Most  of  these  efforts,  however,  proved  fruitless. 

To  begin  with,  a  listing  was  made  of  daily  total  kilowatt 
demand  versus  the  magnitude  of  the  respective  peak  hourly 
load  factor.   It  was  thought  that  a  correlation  might  exist 
which  reflected  higher  peaks  for  those  days  with  higher  over- 
all consumption.   This  was  not  the  case.   Weekdays  with 
relatively  high  electrical  consumption  (average  hourly  demand 
=  10,800  -  11,200  KW)  showed  peaks  ranging  from  1.255  to 
1.391.   Hourly  average  consumption  levels  in  the  range  of 
9600  to  10,000  KW,  however,  had  peak  load  factors  of  1.257 
to  1.4  37.   The  same  type  of  variation  was  present  in  the 
weekend  data.   A  listing  was  also  made  of  average  temperature 
versus  peak  hourly  load  factor  for  each  day.   Again,  no  corre- 
lation was  possible;  the  data  appeared  randomly  distributed. 

Beginning  with  January  1,  1976  and  continuing  in  sequen- 
tial order  through  each  month  of  the  year,  a  matrix  was 
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prepared  for  both  weekends  and  weekdays  in  which  the  hour 
and  magnitude  of  the  hourly  peak  load  factor  were  noted. 
Although  the  magnitude  of  the  peak  varied  widely  from  one  day 
to  another,  the  time  of  occurrence  showed  a  clear  trend  from 
month  to  month.   The  first  three  months  of  the  year  were 
characterized  by  weekday  peaks  at  4:00  P.M.  and  5:00  P.M. 
almost  exclusively.   The  spring  and  summer  months  demonstrated 
mixed  groupings  of  weekday  peaks,  one  around  11:00  or  12:00 
A.M.  and  another  group  at  3:00  P.M.   Peaks  for  the  latter 
three  months  of  the  year  occurred  primarily  in  the  late  after- 
noons with  some  as  early  as  12:00  A.M.   The  data  suggested  that 
the  time  of  occurrence  of  the  daily  peak  was  influenced  by 
ambient  temperature.   A  scatter  plot  of  peak  hour  versus 
temperature  for  the  sampling  period  substantiated  this  trend. 

As  indicated  by  the  fourteen  months  of  data,  the  dominant 
weekday  profiles  for  temperatures  less  than  60 °F  centered 
around  the  hours  of  12:00  A.M.  and  4:00  P.M.   For  temperatures 
above  60°F  a  3:00  P.M.  peak  and  a  12:00  A.M.  peak  were  identi- 
fiable.  Weekends  were  characterized  by  essentially  three 
repeatable  peak  hours.   For  temperatures  less  than  39 °F,  a 
5:00  and  6:00  P.M.  peak  were  dominant.   Between  40°F  and  59 °F 
an  additional  3:00  P.M.  existed.   Above  60 °F  only  a  3:00  and 
5:00  P.M.  peak  were  recorded. 

Inasmuch  as  the  single  most  differentiable  feature  of 
the  electrical  demand  profiles  was  the  time  of  occurrence  of 
the  peak,  the  study  concentrated  upon  finding  the  most 
characteristic  profile  for  each  hourly  peak  grouping.   This 
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analysis  was  made  particularly  difficult  due  to  the  diversity 
of  peak  magnitudes  within  each  grouping.   Bracketed  within  a 
definable  range,  the  magnitudes  of  both  weekday  and  weekend 
peak  electrical  load  factors  appeared  to  be  a  random  variable. 
In  reality,  they  are  not  random  but  rather  a  function  of  campus 
usage.   That  the  data  could  not  be  easily  quantified,  however, 
so  as  to  provide  a  means  for  correlating  peak  magnitudes 
suggested  that  only  the  most  representative  or  typical  peaks 
be  identified. 

For  the  above  purpose,  the  computer  program  previously 
mentioned  in  connection  with  steam  load  profile  determination 
was  employed.   Using  as  input  groups  of  days  which  each 
exhibited  the  same  time  of  peak  occurrence,  the  program  per- 
formed a  least  squares  regression  analysis.   Excellent  results 
were  obtained  despite  the  relative  fluctuation  in  magnitude 
of  the  peak  load  factors.   This  is  attributable  to  the  fact 
that  the  shape  of  the  daily  electrical  profiles  are  very 
similar.   As  with  the  steam  daily  profiles,  some  adjustments 
to  the  polynomial  approximations  were  required  in  order  that 
the  fitted  curves  reflected  accurately  the  data  trends.   For 
example,  almost  every  weekday  profile  showed  a  decreasing 
hourly  load  factor  from  midnight  until  6:00  A.M.   The  higher 
order  polynomial  approximations,  however,  turned  upward  at 
4:00  and  5:00  A.M.,  thereby  not  providing  a  representative 
model.   Similarly,  the  polynomial  approximations  tended  to 
underestimate  the  magnitudes  of  the  daily  peaks,  to  the  extent 
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that  their  use  could  introduce  an  error  of  up  to  10%  in  the 
value  of  the  peak  electrical  demand. 

4.4.1  Weekday  Results 

Five  daily  profiles  were  derived  as  being 
representative  of  typical  weekdays  (Figures  4.12  through  4.16). 
Two  apply  for  temperatures  less  than  or  equal  to  60 °F  and  three 
for  temperatures  greater  than  60 °F.   Because  of  the  spread  in 
magnitude  of  the  daily  peaks  for  days  with  temperatures 
greater  than  60 °F,  two  distinct  profiles  were  specified  for 
days  with  midafternoon  peaks  (Figures  4.14  and  4.15).   The 
following  is  a  breakdown  of  the  proportion  of  days  in  the 
sampling  period  which  exhibited  each  profile. 

Temperature  *  60 °F 

No.  days  with  profile  #1  -  19  (23%) 
No.  days  with  profile  #2  -  63  (77%) 
Total  number  days  in  sample  -  82 

Temperature  >  60°F 

No.  days  with  profile  #1-8  (13%) 
No.  days  with  profile  #2  -  32  (52%) 
No.  days  with  profile  #3  -  22  (35%) 
Total  number  of  days  in  sample  -  62 

4.4.2  Weekend/Holiday  Results 

Three  profiles  were  determined  for  weekends/ 

holidays  (Figures  4.17  through  4.19).   The  proportion  of  days 

in  the  sampling  period  following  each  profile  is  indicated 

below  as  a  function  of  temperature  bandwidth. 
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Temperature  ^  39°F 

No.  days  with  5:00  P.M.  peaks  -  7  (64%) 
No.  days  with  6:00  P.M.  peaks  -  4  (36%) 
Total  number  days  sampled  -  11 

40 °F  *.  Temperature  £  59°F 

No.  days  with  3:00  P.M.  peaks  -  6  (30%) 

No.  days  with  5:00  P.M.  peaks  -  10  (50%) 

No.  days  with  6:00  P.M.  peaks  -  4  (20%) 
Total  number  days  sampled  -  20 

Temperature  2   60 °F 

No.  days  with  3:00  P.M.  peaks  -  19  (70%) 
No.  days  with  5:00  P.M.  peaks  -  8  (30%) 
Total  number  days  sampled  -  27 


4.4.3   Application  of  Daily  Electrical  Load  Profiles 
In  view  of  the  many  factors  which  influence 
the  adherence  to  any  one  particular  daily  electrical  demand 
profile,  it  is  most  difficult  to  establish  the  likelihood  that 
certain  profiles  are  more  dominant  than  others.   Indeed,  the 
data  suggests  almost  a  random  application  of  the  various 
profiles  for  weekdays  and  weekends.   For  the  purposes  of 
electrical  load  simulation,  therefore,  it  is  proposed  that  a 
system  of  random  profile  selection  be  employed.   Since  the 
only  data  examined  regarding  profile  frequency  was  that  for 
January  1976  to  February  1977,  it  has  been  assumed  that  the 
various  load  patterns  apply  for  any  time  frame  in  the  same 
proportion  as  is  reflected  by  the  data  sampling.   For  example, 
within  a  particular  temperature  range,  say  weekdays  less  than 
60°F,  profiles  one  and  two  (Figures  4.12  and  4.13)  can  be 
assigned  to  any  day  providing  that  the  resulting  distribution 
of  these  profiles  is  23%  and  77%  respectively.   A  method, 
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therefore,  exists  for  apportioning  the  representative  consump- 
tion patterns  over  a  range  of  days,  simulating  anticipated 
profiles  using  the  most  well  defined  "curve  fits"  from 
historical  data. 

4.5  Electrical  Load  Profile  Summary 

Tables  4.1  through  4.3  contain  listings  of  the 
hourly  load  factors  for  each  of  the  electrical  load  profiles. 
They  are  intended  as  a  supplement  to  the  graphical  representa- 
tions.  Table  4.4  provides  a  numerical  listing  of  the  daily 
electrical  consumption/temperature  information  which  is 
reflected  by  Figure  4.1.   The  same  information  for  weekends/ 
holidays  (Figure  4.2)  is  shown  in  Table  4.5. 
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Hour 

of 
Day 

Late  Afternoon 
Peaks 

Early  Afternoon 
Peaks 

1 

.820 

.820 

2 

.770 

.780 

3 

.740 

.750 

4 

.715 

.730 

5 

.705 

.720 

6 

.715 

•  715 

7 

.730 

.725 

8 

.840 

.840 

9 

1.060* 

1.060 

10 

1.180 

1.220 

11 

1.230 

1.270 

12 

1.260 

1.300 

13 

1.250 

1.290 

14 

1.270 

1.280 

15 

1.280 

1.270 

16 

1.300 

1.260 

17 

1.250 

1.210 

18 

1.180 

1.130 

19 

1.105 

1.050 

20 

1.010 

.990 

21 

.970 

.950 

22 

.920 

.920 

23 

.870 

.880 

24 

.830 

.840 

Table  4.1  -  Weekday  Electrical  Profile  Hourly  Load  Factors 
for  Days  with  Temperature  £  6o°F. 
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Hour 

of 
Day 

Extreme 

Mid-Afternoon 

Peaks 

Normal 
Mid-Afternoon 
Peaks 

Morning 
Peaks 

1 

.790 

.810 

.800 

2 

.770 

.770 

.770 

3 

.740 

.750 

.750 

k 

.720 

.740 

.740 

5 

.710 

.730 

.730 

6 

.700 

.720 

.720 

7 

.700 

.750 

.740 

8 

.790 

.850 

.860 

9 

1.010 

1.070 

1.080 

10 

1.190 

I.230 

1.220 

11 

1.240 

1.250 

1.310 

12 

1.260 

1.270 

1.310 

13 

1.280 

1.260 

1.260 

14 

1.340 

1.290 

1.290 

15 

1.360 

1.290 

1.280 

16 

1.320 

1.280 

1.270 

17 

1.250 

1.220 

1.210 

18 

1.150 

1.120 

1.110 

19 

1.050 

1.030 

1.020 

20 

1.000 

.980 

.980 

21 

.980 

.950 

.950 

22 

.940 

.920 

.920 

23 

.880 

.880 

.860 

24 

.830 

.840 

.820 

Table 

4.2  -  Weekday  Elec 

trical  Profile  Hour 

ly  Load  Fa< 

for  Days  with  Temperature  >  60°F. 
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Hour 
of 

1500  Peaks 

1700  Peaks 

1800  Peaks 

Day 

1 

.960 

.970 

.960 

2 

.9^5 

.940 

.925 

3 

.930 

.920 

.900 

4 

.915 

.900 

.880 

5 

.900 

.890 

.865 

6 

.870 

.870 

.850 

7 

.850 

.840 

.840 

8 

.880 

.870 

.850 

9 

.940 

.940 

.906 

10 

.990 

.980 

.940 

11 

1.030 

1.020 

.974 

12 

1.060 

1.050 

1.020 

13 

1.100 

1.070 

1.050 

14 

1.100 

1.090 

1.075 

15 

1.130 

1.100 

1.100 

16 

1.110 

1.105 

1.120 

17 

1.100 

1.140 

1.150 

18 

1.070 

1.105 

1.170 

19 

1.030 

1.090 

1.140 

20 

1.040 

1.060 

1.120 

21 

1.050 

1.040 

1.098 

22 

1.030 

1.030 

1.064 

23 

1.000 

1.000 

1.024 

24 

.970 

.980 

.979 

Table  4.3  -  Weekend/Holiday  Electrical  Profile  Hourly 
Load  Factors. 
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V     LOAD  GROWTH  AT  MIT 

Consideration  of  total  energy  system  design  must  neces- 
sarily take  into  account  future  campus  energy  needs.   Plant 
sizing  requires  a  realistic  estimate  of  the  projected  loads 
associated  with  Institute  growth  over  some  finite  time  span. 
An  assessment  of  this  sort  can  be  strongly  influenced  by 
measures  taken  to  limit  energy  usage,  above  and  beyond  those 
associated  with  innovative  building  design  and  construction. 
For  this  reason  it  is  worthwhile  to  examine  methods  which  have 
been  employed  during  recent  years  at  MIT  in  connection  with 
energy  conservation  and  define  the  trend  for  the  future. 

5 . 1   Overall  Campus  Energy  Conservation  Measures [ 12 ] 

As  early  as  1969 ,  MIT  began  actively  seeking  ways 
to  use  energy  more  efficiently.   In  that  year  a  power  factor 
correction  program  was  instituted.   Prior  to  this  time  a 
relatively  low  power  factor  for  campus  electrical  usage  was 
in  part  responsible  for  higher  billing  charges  from  Cambridge 
Electric.   At  a  cost  of  approximately  $60,000,  a  system  of 
capacitors  was  installed  for  the  purpose  of  increasing  the 
power  factor.   This  program  was  entirely  paid  for  prior  to 
more  serious  energy  conservation  measures  were  undertaken  in 
the  early  1970' s. 

In  1972,  MIT  became  aware  that  the  incremental  cost  of 
chilled  water  from  the  Central  Plant  (fuel  fired  only)  was 
considerably  lower  than  the  cost  of  electric  air  conditioning 
in  a  number  of  campus  buildings.   It  was,  therefore,  decided 
to  connect  these  buildings  to  the  Central  Plant  chilled  water 
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mains.   The  cost  of  connection  was  recovered  in  a  two  year 
payoff  program.   Conservation  activities  in  limiting  air 
conditioning  loads  across  the  Central  Utilities  Plant  service 
area  were  successful  in  reducing  the  chilled  water  demand  by 
approximately  2000  tons  from  1973  to  1974. 

Sparked  by  the  Arab  oil  embargo  of  1973,  emergency 
measures  were  taken  to  reduce  campus  fuel  consumption.   The 
temperature  of  occupied  areas  in  MIT  buildings  was  set  back 
to  68°F  during  the  heating  season  while  that  in  unoccupied 
areas  was  reduced  to  a  maximum  of  50 °F.   All  air  conditioning 
except  that  used  for  computer  and  experimental  work  was  cur- 
tailed during  the  winter  months.   Hot  water  supply  temperature 
was  lowered  from  140 °F  to  100 °F  in  lavatories  while  that  for 
dormitories  was  reset  to  120 °F.   Such  conservation  efforts 
continue  today  in  the  face  of  rising  energy  costs. 

In  addition  to  the  above,  significant  lighting  reductions 
have  been  achieved.   Corridor  lighting  has  been  set  to  a  level 
of  five  foot  candles  (fc)  for  normal  use.   Office,  laboratory 
and  classroom  lighting  levels  have  been  specified  to  be  no 
greater  than  50  to  70  fc.   All  decorative  and  outdoor  archi- 
tectural lighting  has  been  eliminated.   A  program  for  convert- 
ing incandescent  to  fluorescent  (or  other  high-efficiency) 
illumination  has  been  instituted. 

5.2  Load  Management  at  MIT 

Control  of  campus  energy  usage  presently  depends  upon 
the  efforts  taken  by  both  man  and  machine.   At  the  Central 
Utilities  Plant,  watch  personnel  frequently  switch  from  steam 


-124- 

to  electric  driven  auxiliaries  and  back  again  so  to  achieve 
the  most  efficient  equipment  mix  for  the  prevailing  campus 
demands.   Time  clocks  in  many  buildings  provide  for  scheduled 
start  and  stop  of  heating,  ventilation  and  air  conditioning 
machines.   Peak  shaving  with  a  1100  KW  generator  is  undertaken 
on  a  regular  basis  during  the  higher  demand  period  of  each 
weekday.   Of  more  significance,  however,  is  the  role  played 
by  computers  [  3  ]  in  the  reduction  of  power  consumption  at  MIT. 

In  reaction  to  the  energy  crisis  of  1974,  MIT  made  a 
preliminary  study  of  the  cost  benefits  which  could  be  achieved 
using  a  central  computer  controller  for  power  management 
purposes.   Indications  were  that  although  such  a  system  was 
relatively  expensive  to  purchase,  it  showed  a  high  rate  of 
return.   As  a  result,  an  IBM  System/7  sensor-based  computer 
was  installed  and  brought  on  line  in  December  of  1974. 
Connected  to  eight  buildings  with  unusually  large  energy 
consumption  levels,  the  computer  provided  demand  control 
through  equipment  cycling  as  well  as  peak  load  clipping. 
During  its  period  of  use  the  IBM  System/7  computer  more  than 
paid  for  itself  as  savings  rates  of  ten  percent  were  realized 
in  the  buildings  monitored.   Moreover,  it  served  a  proving 
ground  function  inasmuch  as  it  provided  the  Department  of 
Physical  Plant  with  sufficient  confidence  in  the  concept  of 
automated  central  controls  that  a  more  sophisticated  and  wide- 
spread computer-based  management  system  was  ultimately  adopted, 

Known  as  the  Facilities  Management  Systems  (FMS) ,  the 
present  means  of  effecting  load  control  is  through  a  central 


-125- 

computer  station  utilizing  dual  PDP-11/40  processors.   The 
system  communicates  with  remote  stations  located  in  thirty- 
four  campus  buildings.   Within  each  building  the  remote 
station,  in  turn,  communicates  with  various  sensing  and 
actuating  elements  and  reports  back  to  the  central  computer 
station.   Optimized  equipment  start/stop  schedules  are  moni- 
tored by  the  system,  and  HVAC  units  are  managed  so  as  to  pro- 
vide a  desirable  level  of  humidity  and  temperature  control. 
Dampers,  humidifiers,  cooling  coils  and  heating  coils  are 
directly  scheduled  by  the  FMS.   Every  fifteen  minutes  the 
current  position  of  each  HVAC  system  valve  or  actuator  is 
monitored.   Variations  in  building  loads  are  adjusted  for 
automatically  so  that  heating  or  air  conditioning  distribution 
loops  just  satisfy  demand,  thereby  conserving  energy  and 
dollars. 

Electrical  consumption  is  managed  by  cycling  loads  and 
limiting  the  peak  demand  during  a  day.   Cycling  of  building 
loads  is  accomplished  using  a  computer-controlled  schedule 
for  the  starting  and  stopping  of  "sheddable"  (non-essential) 
electrical  loads.   Equipment  so  managed  is  cycled  ON  and  OFF 
many  times  per  day,  with  a  variable  OFF  interval  between  five 
and  sixty  minutes  per  hour. 

Limitation  of  peak  electrical  demands  is  a  two  phase 
operation.   As  demand  increases,  the  central  computer  station 
generates  an  advisory  message  to  the  Central  Utility  Plant  to 
bring  the  1100  KW  emergency  generator  on-line.   As  a  supple- 
ment to  the  normal  electrical  supply  from  Cambridge  Electric, 
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use  of  this  generator  serves  to  reduce  the  peak  monthly 
kilowatt  load  on  which  demand  charges  are  computed.   In 
addition  to  the  above,  load  shedding  is  undertaken  as  a  means 
of  reducing  peak  demands.   The  entire  campus  electrical  load 
is  constantly  measured.   Based  on  the  load  values  for  four 
consecutive  three  minute  intervals,  a  kilowatt  forecast  is 
made  for  the  next  thirty  minutes  using  a  linear  regression 
technique.   This  forecast  is  compared  to  a  preset  demand 
limit  target.   If  it  appears  that  the  target  will  be  exceeded, 
shedding  of  non- crucial  loads  is  implemented.   When  the 
forecast  no  longer  predicts  that  the  demand  limit  will  be 
exceeded,  load  shedding  ceases.   The  hierarchy  of  priority 
for  load  reduction  with  the  FMS  is: 

(a)  Load  Shedding  Procedures 

(b)  HVAC  Optimization 

(c)  Load  Cycling 

(d)  Optimized  Start/Stop 

At  present,  MIT's  FMS  monitors  buildings  which  total 
about  four  million  square  feet  and  account  for  75%  of  total 
energy  consumption.   It  controls  approximately  2500  points. 
The  communications  trunks  have  the  capacity  for  several  times 
the  traffic  they  now  are  carrying.   Although  the  computer 
subsystem  can  be  expanded,  it  can  accommodate  two  to  three 
times  the  number  of  transactions  it  now  processes.   Aside  from 
the  requirement  of  adding  remote  stations  to  place  additional 
buildings  under  FMS  control,  there  is  no  practical  restriction 
on  growth  of  the  system. 
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Through  the  central  control  of  some  2500  different  valves 
or  actuators  in  thirty-four  different  buildings,  a  campus-wide 
integrated  systems  approach  has  been  realized  with  FMS.   As 
more  buildings  join  the  management  grid,  the  control  of  campus 
energy  usage  will  become  more  exact.   Even  now,  in  the  brief 
time  FMS  has  been  operational,  the  "single  system"  management 
approach  has  yielded  efficiencies  which  were  unobtainable 
before. 

Present  estimates  are  that  once  FMS  is  extended  to  every 
campus  building,  load  reductions  (steam  and  electric)  on  the 
order  of  12%  will  be  realized.   This  figure  is  relative  to 
current  levels  of  energy  consumption,  prior  to  implementation 
of  FMS. 

5. 3   Long  Range  Building  Plans  at  MIT 

Available  through  the  Office  of  Planning  is  tentative 
information  on  both  the  type  and  square  footage  of  buildings 
MIT  is  likely  to  fund  for  construction  through  the  year  2000. 
While  most  of  the  plans  are  for  new  construction,  renovation 
of  existing  buildings  accounts  for  approximately  10%  of  the 
development  space. 

For  presentation  herein,  buildings  have  been  grouped 
according  to  type.   Five  classifications  are  identified: 

(a)  Classroom/Faculty  Office  Building 

(b)  Classroom/Laboratory/Workshop  Building 

(c)  Administrative  Offices 

(d)  Residential  Building 
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(e)   Athletic  Building 

Table  5.1  summarizes,  in  five-year  increments,  the  contribution 
of  each  building  group  to  the  total  area  which  is  scheduled. 

Of  importance  to  this  study  is  the  fact  that  every 
building  addition  which  is  planned  up  to  each  of  the  five 
year  milestones  represents  new  steam  and  electrical  demands  on 
a  potential  total  energy  system.   Although  some  renovation  is 
included  in  the  plans,  the  buildings  so  designated  are  not 
presently  on  the  campus  grid.   That  is,  they  are  billed 
separately  for  the  electricity  they  consume  and  do  not  form 
any  part  of  what  is  currently  recognized  as  MIT's  normal 
electrical  load. 

Information  on  campus  growth  is  of  little  use  without  a 
meaningful  breakdown  of  the  steam  and  electrical  loads 
associated  with  the  respective  types  of  buildings  MIT  plans  to 
construct. 

5. 4   Load  Estimation  for  Future  Campus  Construction 

For  the  purposes  of  total  energy  system  design  and 
selection  it  is  desirable  that  an  estimation  be  made  as  to  the 
probable  extent  of  the  load  increase  attributable  to  the  long 
range  building  program  at  MIT.   As  this  impacts  significantly 
upon  equipment  selection  and  sizing,  such  an  estimation  must 
not  assume  that  the  past  is  likely  to  be  repeated  where  build- 
ing energy  consumption  levels  are  concerned. 
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1980      1985      1990      1995      2000 

Classroom/Faculty 

Office  Building 

Sq.  ft.      —     250,000    150,000    600,000    100,000 
to  Date     "~     250,000   400,000   1,000,000  1,100,000 

Classroom/Laboratory/ 
Workshop  Building 

Sq.  ft.   305.000    285,000    110,000    160,000      — 
to^Date   305.000   590,000   700,000   860,000   860,000 

Administrative 

Offices 

Sq.    ft.       115,000        135.000  —  —  _ 

to^Date      n5.000        250,000        250,000        250,000        250,000 

Residential 

Sq.  ft.     —     120,000   295.000   230,000   190,000 

—    120,000   415,000   645,000   835,000 


Total 
to  Date 


Athletic 
Building 

Sq.  ft.      —     225,000      —        —        — 
to^Date     ""     225,000   225,000   225,000   225,000 

Total  Sq.  ft. 

Added  During 

Previous  Five 

Years         420,000   1,015,000   555,000   990,000   290,000 

Total 

to  Date       420,000   1,435,000  1,990,000  2,980,000  3.270,000 


Table  5.1  -  MIT  Long  Range  Building  Plans  by  Building  Type 
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5.4.1   Intensity  of  Energy  Usage;   1960  -  1976 

If  the  trend  in  new-building  energy  consump- 
tion during  the  past  fifteen  years  provides  any  indications 
for  the  future,  one  of  them  is  that  new  building  design  must 
undergo  radical  changes  immediately.   A  1974  study  by  the  MIT 
Environmental  Engineer  showed  that  academic  buildings  con- 
structed between  1960  and  1970  have  electrical  use  levels 

2 

ranging  from  10.5  to  54  KWh/ft  -year  with  a  median  of 

2 
35  KWh/ft  -year.   In  contrast,  the  original  buildings  of  MIT's 

2 

main  group  consume  electricity  at  a  rate  of  10.5  KWh/ft  -year. 

The  disparity  in  the  two  consumption  figures  is  due  to  several 
factors  [  l  ] ,  among  which  are  the  large  air  circulation 
requirements  for  research  hoods  and  increased  modular  lighting 
levels  of  newer  buildings.   Trends  in  dormitory  construction 
have  done  nothing  but  aggravate  an  already  worsening  situation. 
High-rise  building  design  and  the  replacement  of  commons  dining 
with  individual  kitchens  have  doubled  the  intensity  of 
electrical  energy  use  over  the  traditional,  older  dormitory 
style. 

Consumption  of  heat  energy  has  followed  patterns  similar 
to  those  mentioned  above  in  buildings  of  more  recent  construc- 
tion.  Actual  data  shows  that  the  design  of  newer  buildings  is 
such  that  the  yearly  heating  load  is  approximately  twice  that 
of  the  older  buildings. 

That  any  future  campus  buildings  will  incorporate  energy 
efficient  designs/systems  is  a  foregone  conclusion.   The  only 
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unknown  is  the  degree  of  load  reduction  which  can  be  achieved 
by  good  design  alone. 

5.4.2   Examination  of  Present  Usage  Data 

In  determining  a  plausible  estimate  for 
future  steam  and  electrical  usage  levels,  consumption  figures 
for  existing  MIT  buildings  were  reviewed.   Buildings  were 
grouped  according  to  the  five  types  which  were  addressed  in 
Section  5.3.   For  each  group  of  buildings,  data  showing  the 
present  intensity  of  steam  and  electrical  energy  use  was 
listed.   Based  on  these  figures  estimates  were  made  as  to 
possible  "design"  consumption  levels  which  might  be  achieved 
in  the  future.   The  information  is  summarized  in  Table  5.2. 

The  buildings  chosen  for  inclusion  in  Table  5.2  comprise 
a  mix  of  newer  and  older  construction.   Therefore,  their 
energy  consumption  levels  show  considerable  variance.   It  can 
be  argued  that  if  buildings  of  low  energy  consumption  were 
built  in  the  past,  they  can  be  built  again.   Indeed,  in  the 
face  of  growing  concern  over  energy  costs  this  must  occur. 
Reliable  figures  are  not  available,  however,  on  which  to  base 
definitive  intensity-of-use  estimates  for  the  future  construc- 
tion of  each  building  type.   The  numbers  presented  in  Table 
5.2  reflect  the  author's  belief  that  the  most  realistic 
values  lie  between  the  upper  and  lower  extremes,  with  the 
lower  extreme  being  more  heavily  favored.   In  any  case,  the 
new  building  consumption  levels  are  only  estimates  to  be  used 
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Building  Type  Intensity  of  Use 

Electrical        Steam 

(KWh/ft2-yr)  (lbs/ft^-yr) 

72-73   75-76  72-73   75-76 

Classroom/Faculty 

Office  Bulljlnjc 

Sloan  (E52)  . 12.785   10.222  107.0   65.25 

Pairchild( 36/38)  24.008  18.580  97.6   67.4 

Compton  (26)  16.764  10.911  170.4   60.9 

Space  Research  (37)  ...  16.584  14.726  102.3   61.4 

Future  Design  Estimate       12.0  62.0 

Classroom/Laboratory/ 

Workshop  Building 

Bush  113) 42.214  25.567  228.7  137.1 

Dorrance  (16)  22.083  20.206  70. 8   61.4 

Whitaker  (56)  54.377  28.120  182.8  109.6 

Future  Design  Estimate       23. 0  80. 0 

Administrative 

Offices 

Ford  (E18)  21.488  16. 450  142.2  135.6 

Maclaurin  (10)  IO.613   6.389  102.4   61.4 

Future  Design  Estimate       10.0  80.0 

Residential 

Eastgate  (E55)  9.720   8.382  109.4   84.6 

Baker  (W7)  4.672   4.040  143.3  126.1 

Tang  (W84)  7.904   7.36O  124.9   85. 0 

McCormick  (W4)  10.479   8.253  101. 3   64.4 

MacGregor  (W6l)  9.506   7.139  96.3   81. 5 

Burton  (W51)  6.330   5.^54  78.1   75. ^ 

Future  Design  Estimate       7*0  85. 0 

Athletic 
Building 
Hfepont  (W32)  7.53*   ^536  87.4   53.8 

Future  Design  Estimate       5*5  65*0 


Table  5.2  -  Usage  Intensity  Information  for  Selected  Types 
of  MIT  Buildings 
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in  the  relative  sizing  of  plant  capacity  to  meet  campus 
steam  and  electrical  needs  in  the  future. 

5.5  Projected  Institute  Electrical  &  Steam  Load  Growth 
Using  the  information  from  Tables  5.1  and  5.2,  pro- 
jections of  campus  load  growth  were  made  for  each  five  year 
time  frame  commencing  in  1980. 

The  average  kilowatt  increase  resulting  from  new  build- 
ing construction  was  computed  as, 


(square  footage] 
added 


estimated  electrical  usage  intensity, 
KWh/ft2-year 


Similarly,  the  average  steam  demand  increase  was  computed  as 


[square  footage] 
added 


estimated  steam  usage  intensity , 
2 
lbs/ft  -year 


By  themselves,  average  hourly  consumption  figures  are  diffi- 
cult to  interpret.   A  meaningful  reference  is  provided  by  1976 
consumption  information. 

For  1976  the  average  hourly  electrical  load  was  9737 
kilowatts  (total  kilowatt  hours  used/8760  hours) .   Average 
hourly  steam  consumption  was  89,419  lb/hr.   The  peak  electri- 
cal load  in  1976  was  15,240  KW,  representing  a  demand  1.56 
times  as  great  as  the  hourly  average  electrical  load.   The 
peak  steam  load  in  1976  was  217,000  lbs,  reflecting  a  heating 
demand  2.43  times  greater  than  the  hourly  average.   If  these 
same  factors  are  applied  to  each  five  year  step  increase  in 
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campus  load,  an  appreciation  for  the  required  plant  capacity 
to  meet  peak  demands  can  be  gained.  Tables  5.3  and  5.4  pro- 
vide this  summary. 

It  can  be  seen  that  in  order  to  just  meet  the  projected 
demands  in  1985,  a  total  energy  system  which  can  accommodate 
a  19,000  KW  load  and  a  250,000  lb/hr  steam  load  must  be 
installed.   By  1990  these  figures  increase  to  20,100  KW  and 
260,000  lb/hr  steam.   Allowing  for  a  design  margin  of  20%  (to 
absorb  unforeseen  load  growth  past  1990)  the  above  demands 
dictate  that  an  electrical  generation  capacity  of  24,120  KW 
be  installed  with  a  combined  boiler  capacity  of  312,000  lb/hr. 

The  above  considerations  ignore  the  role  FMS  plays  in 
reducing  the  daily  peak  loads.   Using  the  most  recent  infor- 
mation available  on  the  system  performance,  a  3%  peak  load 
reduction  can  be  assumed  to  exist.   It  is  recalled  from 
Section  5.2  that  load  reductions  of  up  to  12%  a  year  are 
anticipated  once  FMS  is  fully  integrated  into  all  campus  build- 
ings.  Allowing  for  this  magnitude  of  load  reduction  lowers 
the  average  hourly  steam  and  electrical  demands  projected  for 
1990  to  115,150  lb/hr  and  13,912  KW  respectively. 

Table  5.5  summarizes  the  information  pertinent  to  sizing 
of  a  total  energy  system  to  accommodate  demands  in  1990.   It 
reflects  calculations  for  a  design  margin  of  20%  as  well  as 
load  reductions  attributable  to  FMS.   It  shows  that  an  elec- 
trical generating  capacity  of  21,050  KW  is  needed  and  a 
boiler  capacity  of  271,425  lb/hr.   These,  then,  can  be  used 
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Projected  Campus  ^   . 

Steam  Demand    (lb/yr)    93^.012,138 

12#  Overall  Load  ._  ,_     . 

Reduction   (FMS)      "  U2.08l.457 

20£  Design  Margin     +  186,802.428 

Net   ....    1,008,733.109 

Hourly  Average 

Demand        115.152 

Projected  Peak  279,820 

3^  Peak  Reduction  -  8,395 

Net  271,^25 


^uct^Uf -13.5^0.620 

20#  Design  Margin  +  22,567,700 

Net 121,865,580 

Hourly  Average  n  oi? 

Demand        0»V" 

Projected  Peak  21,700 

3*  Peak  Reduction -  650 

Net  21,050 


Table  5-5  -  Projected  Peak  Demands  for  the  Sizing  of 
Plant  Equipment  to  Satisfy  1990  Loads. 
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as  design  figures  for  the  sizing  of  equipment  in  proposed 
total  energy  system  schemes. 

It  has  been  assumed  that  steam  heating  and  air  condition- 
ing will  be  provided  from  the  Central  Utility  Plant  for  all 
new  building  construction.   This  is  not  a  binding  requirement 
It  may  not  be  feasible  to  design  every  new  building  addition 
so  that  it  is  a  "natural"  load  extension  of  a  centrally 
located  total  energy  plant.   In  the  place  of  chilled  water 
from  a  campus  wide  loop,  therefore,  design  engineers  at  some 
future  time  may  elect  to  place  electric  driven  compressors 
in  a  new  building  to  avoid  overloading  the  Central  Chiller 
Plant.   Such  decisions  will  be  made  on  a  case  basis.   That 
some  latitude  must  exist  in  the  forecasting  of  the  specific 
mix  of  load  additions  is  a  necessity  if  an  optimum  thermal 
to  electric  load  ratio  is  to  be  achieved. 
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VI    REPRESENTATIVE  YEAR  MODEL 

As  a  result  of  the  work  completed  in  Chapters  III  and  IV, 
a  methodology  exists  by  which  to  model  both  steam  and  electri- 
cal demands  at  MIT.   The  motivation  for  this  development  has 
been  the  desire  to  simulate  campus  energy  requirements  in 
various  total  energy  system  schemes  for  the  purpose  of  deter- 
mining the  relative  cost  advantages  of  each.   Of  primary 
interest  in  this  investigation  is  the  efficiency  with  which 
a  particular  equipment  configuration  satisfies  both  steam  and 
electrical  demands. 

Where  a  severe  mismatch  between  hourly  thermal  and 
electric  loads  exists,  certain  plant  designs  appear  much  more 
attractive  than  others.   Conversely,  for  demand  patterns  whose 
shapes  track  closely,  the  economics  of  selection  dictate  that 
still  other  designs  are  preferred.   In  the  case  of  MIT, 
demand  profiles  conform  to  no  one  pattern  but,  instead, 
fluctuate  considerably  during  the  course  of  a  year.   For  this 
reason,  the  comparison  of  different  plant  designs  must  be 
based  on  data  which  is  representative  of  the  entire  spectrum 
of  thermal  to  electric  load  ratios. 

Chapters  III  and  IV  demonstrated  the  importance  of  ambient 
temperature  as  a  parameter  for  predicting  both  daily  total 
steam  and  electrical  requirements  at  MIT.   To  be  certain,  the 
reliability  of  each  load  model  is  heavily  dependent  upon  the 
availability  of  valid  yearly  temperature  information  for  the 
Boston  area.   Had  1976  been  a  truly  representative  year  as 
far  as  temperature  is  concerned,  the  load  data  for  that 
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period  could  be  used  directly  in  a  computer  simulation  as 
being  "typical"  of  any  average  year.   Such  was  not  the  case. 
The  first  three  months  of  winter  were  between  three  and  four 
degrees  cooler  than  usual  while,  the  late  winter  and  spring 
months  were  approximately  five  degrees  warmer  than  usual.   Con- 
sequently, it  was  required  to  develop  a  "model  year"  in  order 
to  form  a  basis  for  load  inputs  into  a  computer  simulation 
program. 

6.1   Data  Collection 

Available  at  the  Boston  Weather  Bureau  is  historical 
temperature  data  for  every  day  of  the  year.   Table  6.1  sum- 
marizes some  of  this  information.   For  each  day  a  seventeen 
year  smeared  average  of  the  daily  mean  observed  temperature 
is  listed.   The  data  is  presented  relative  to  a  12.5  mph 
prevailing  wind  in  the  Boston  area.   While  a  considerable 
temperature  distribution  can  be  observed  from  winter  to 
summer,  the  large  number  of  days  in  January  and  February  with 
mean  temperatures  at  30°F  is  misleading.   The  direct  result  of 
smoothing  out  temperature  extremes,  this  grouping  renders  a 
direct  usage  of  Table  6.1  impractical. 

Table  6.2  incorporates  the  same  data  as  Table  6.1  but  in 
different  form.   As  an  aid  in  determining  the  evenness  of 
temperature  distribution  over  a  year,  it  focuses  attention  on 
the  upper  temperature  extreme.   Taken  literally,  it  suggests 
that  a  normal  year  in  Boston  has  no  days  with  an  average 
temperature  greater  than  76 °F.  Again  the  result  of  a  seventeen 
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Month   JAN  FEB  MAR  APR  MAY  JUN  JIJL  AUG  SEP  QHT  Nov  http 


Day 
1 

2 

3 

4 

5 
6 

7 

6 

9 

10 

11 

12 

13 

14 

15 

16 

17 

16 

19 

20 

21 

22 

23 
24 
25 
26 
27 
28 
29 
30 
2L 


30   30   33  42  54  65  71   75  68  62  49  40 


31 


32 


37 


41 


36  46 

I 
47 


38 
J   49 
39 

50 

I 

51 

40 


52 


42  5#f 

i 


34  43 

1 

44  55 

I   J 

45  56 


57 


48  I 
59 


60 


61 


53  64 


11 


*    1 

66  72 


67 


62  70 


63 


I   68  74 

58 


76 


71  75 


73  74 


73 


69  75   J 
71 


67 


61 
60 


59 


58 


66 


7.0  64 


69  63 


68 

1  * 


53 


52 


39 

38 

I    I 

48  37 

i 


47 


57     46 


56 


43 


51  42 


50  41 


40 


36 


35 

I 
34 


I   ^5  33 
55 
54 

44  32 


31 


30 


Table  6.1  -  Historical  Daily  Average  Temperatures  (*P) 
for  Boston  Area. 
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-ature  (°f; 

)   #  of  Days 

Temperature  (°F) 

//  of  Days 

30 

59 

54 

6 

31 

8 

55 

3 

32 

6 

56 

6 

33 

4 

57 

6 

34 

6 

58 

4 

35 

6 

59 

6 

36 

6 

60 

6 

37 

6 

61 

3 

38 

4 

62 

6 

39 

6 

63 

5 

40 

6 

64 

7 

41 

5 

65 

9 

42 

5 

66 

10 

43 

5 

67 

11 

44 

6 

68 

11 

45 

6 

69 

9 

46 

6 

70 

11 

47 

6 

71 

7 

48 

6 

72 

7 

49 

7 

73 

9 

50 

6 

74 

8 

51 

6 

75 

19 

52 

5 

76 

3 

53 

5 

Table  6.2  -  Yearly  Listing  of  Number  of  Days  Versus  Average 
Temperature  from  Historical  Weather  Data  for 
Boston. 
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year  averaging,  the  large  grouping  of  days  at  75 °F  is  hardly 
representative  of  any  one  particular  year  in  Boston. 

6. 2  Temperature/Seasonal  Model 

With  the  exception  of  the  temperature  extremes 
mentioned  above,  the  remaining  spread  of  temperatures  is  quite 
evenly  distributed  over  the  year.   For  use  in  a  temperature 
model  it  is  assumed  that  this  median  range  is  representative 
of  a  typical  year  in  Boston.   A  review  of  temperature  records 
for  the  past  ten  years  supports  this  assumption.   A  method 
is  needed,  however,  to  apportion  the  days  with  extreme  tempera- 
tures (both  high  and  low)  over  a  more  realistic  bound. 

6.2.1   Refinement  of  Temperature  Distribution 

In  order  to  gain  insight  into  the  number  of 
days  with  temperatures  less  than  30 °F,  data  for  each  of  the 
last  ten  years  was  examined.   Listings  were  made  of  daily 
average  temperature  and  the  number  of  days  within  each  year 
at  that  temperature.   For  the  ten  year  period,  summations  were 
made  of  the  total  number  of  days  at  each  temperature  less  than 
30 °F  .   It  was  found  that  364  days  exhibited  average  tempera- 
tures in  that  category.   Using  this  figure  as  a  base,  percent- 
age distributions  of  days  at  each  temperature  were  computed 
over  the  entire  range.   The  percentages  were  converted  to 
number  of  days  through  multiplication  by  59  (the  number  of 
days  at  30°F  from  Table  6.2).   The  end  result  was  a  spread  of 
temperatures  extending  to  8°F.   A  more  valid  apportionment 
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than  simply  59  days  at  30°F,  this  refined  temperature  distri- 
bution was  taken  as  being  representative  of  a  typical  year  in 
Boston. 

Essentially  the  same  procedures  were  employed  for  tempera- 
tures at  the  upper  extreme.   Data  was  examined  over  a  ten  year 
period  for  temperatures  greater  than  65°F  and  percentage 
distributions  were  again  computed.   Multiplication  by  107  (the 
number  of  days  greater  than  65°F  from  Table  6.2)  yielded  the 
number  of  days  at  each  temperature. 

Table  6.3  shows  the  final  result  of  using  a  ten  year 
cross  section  for  both  temperature  extremes.   This  distri- 
bution of  temperature  was  the  basis  for  the  model  year. 

6.2.2   Seasonal  Temperature  Breakdown 

Because  the  daily  steam  load  profiles  are 
characteristic  of  specific  seasons,  it  was  necessary  to  deter- 
mine an  approximate  dividing  line  for  temperature  within  each 
season.   For  this  purpose  the  steam  data  logs  at  the  Central 
Utility  Plant  were  reviewed. 

Daily  profiles  were  examined  for  several  weeks  within 
each  of  the  winter,  spring,  summer  and  fall  periods.   The  first 
decidedly  spring  profile  was  found  to  coincide  roughly  with 
the  start  up  of  the  Central  Utility  Plant's  chiller  system. 
In  1976  this  took  place  during  mid-April.   Conversations  with 
Mr.  George  Reid,  Assistant  Chief  of  the  Plant,  revealed  that 
no  set  date  exists  for  beginning  warm  weather  operation  of 
the  system.   He  suggested,  however,  that  mid-April  was  typical. 
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Temp    (°F)      //  Days       Temp    (°F)      #  Days       Temp    ("F)      # 


Days 


8 

35 

6 

61 

3 

10 

36 

6 

62 

6 

11 

37 

6 

63 

5 

12 

38 

4 

64 

7 

13 

39 

6 

65 

6 

14 

40 

6 

66 

5 

15 

41 

5 

67 

6 

16 

42 

5 

68 

6 

17 

2 

43 

5 

69 

6 

18 

3 

44 

6 

70 

8 

19 

2 

45 

6 

71 

9 

20 

3 

46 

6 

72 

6 

21 

4 

47 

6 

73 

5 

22 

2 

48 

6 

74 

6 

23 

4 

49 

7 

75 

5 

24 

3 

50 

6 

76 

6 

25 

4 

51 

6 

77 

6 

26 

4 

52 

5 

78 

4 

27 

4 

53 

5 

79 

6 

28 

4 

54 

6 

80 

4 

29 

5 

55 

3 

81 

3 

30 

7 

56 

6 

82 

4 

31 

8 

57 

6 

83 

3 

32 

6 

58 

4 

84 

1 

33 

4 

59 

6 

85 

2 

34 

6 

60 

6 

Table  6.3  -  Final  Model  Year  Temperature  Distribution 
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As  moderate  to  heavy  use  of  the  chiller  system  accompanies  the 
hot  summer  weather,  so  too  does  the  daily  steam  load  profile 
change.   It  was  found  that  an  approximate  temperature,  both 
in  late  spring  and  early  fall,  which  bracketed  the  heavier 
usage  period  was  65 °F.   No  firm  cut  off  date  exists  for 
shut  down  of  the  chiller  system.   For  1976,  operation  continued 
into  November.   The  experience  of  recent  years,  however, 
suggests  that  late  October  is  more  typical. 

Based  on  the  above,  four  temperature  bandwidths  were 
constructed  which  reflect  an  approximate  breakdown  by  season. 
Using  Table  6.1  as  a  guide  in  assigning  dates,  the  following 
time  frames  were  chosen  as  being  representative  of  "model 
seasons" : 

Winter:   November  1  -  April  14  (165  days) 

(49°F  48°F) 

Spring:   April  15  -  May  31  (  47  days) 

(48°F       65°F) 

Summer:   June  1  -  September  15  (107  days) 

(65°F  66°F) 

Fall:   September  16  -  October  31  (  46  days) 

(65°F  50°F) 

For  modeling  purposes,  assignment  of  a  specific  temperature 
to  days  within  any  one  season  was  made  in  accordance  with 
Table  6.3  and  the  seasonal  listings  above. 

6. 3  Weekday/Weekend  Temperature  Assignment 

The  development  of  a  "typical"  year  implies  that 
2/7  of  the  days  are  weekends  while  5/7  are  weekdays.   The 


•1^7- 


consideration  of  holidays,  however,  upsets  this  balance 
slightly.   Twelve  holidays  are  recognized  for  MIT  employees 
during  a  normal  year,  distinguished  from  that  group  of  student 
holidays  which  do  not  include  the  entire  MIT  community.   When 
added  to  the  104  weekend  days  this  yields  a  total  of  116  days 
which  exhibit  weekend/holiday  demand  characteristics.   249 
weekdays  remain. 

The  assignment  of  specific  temperatures  to  weekdays  and 
weekends/holidays  was  made  as  follows: 

(a)  For  every  season  a  numerical  listing  was  made  of 
each  day  and  the  temperature  ascribed  to  it  (from 
Section  6.2.2)  . 

(b)  Beginning  with  the  first  day  of  winter  five  weekdays 
were  specified,  followed  by  two  weekends.   This 
procedure  was  repeated  throughout  the  year  until  all 
the  days  were  assigned. 

(c)  The  above  listing  was  modified  to  reflect  the 
proper  number  of  holidays  within  each  season.   That 
is,  where  a  season  had  too  few  holidays,  an  appro- 
priate number  of  weekdays  would  be  deleted. 

The  result  of  the  above  assignment  procedure  was  a  grouping 

ot  weekdays  and  weekends/holidays  by  season,  each  with  a 

specific  temperature  ascribed  to  it  (Tables  6.4  -  6.7).   The 

next  step  in  modeling  the  representative  year  was  the  detailing 

of  total  daily  steam  and  electrical  loads  to  each  day. 

6.4   Daily  Demand  Assignment 

A  straightforward  application  of  equations  3.1  and 
3.2  yields  the  daily  total  steam  load  for  each  weekday  and 
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Temperature  (°F) 


Weekdays 


Weekends/Holidays 


8  1 

10  1 

11  1 

12  1 

13  1 

14 o 

15  o 

16  l 

17  2 

18  2 

19  1 

20  2 

21  3 

22  1 

23  3 

24 2 

25  3 

26  3 

27  3 

28  4 

29  3 

30  4 

31  5 

32  4 

33  2 

34  4 

35  4 

36  4 

37  4 

38  3 

39  4 

40  4 

41  3 

42  4 

8  :::::::::::::::  2 

45  4 

46  4 

47  4 

48  2 

49  3 


Total 


111 


J* 


Table  6.4  -  Weekday  and  Weekend  Temperature  Distribution  for 
Model  Winter. 
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Temperature  (°F) 

48  

49  

50  

51  

52  

53  

54  

55  

56  

57  

58  

59  

60  

61  

62  

63  

64 

65  


Weekdays 


Weekends/Holidays 


Total 


32 


Table  6.5  -  Weekday  and  Weekend  Temperature  Distribution  for 
Model  Spring. 
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Temperature  (°F)      Weekdays      Weekends/Holiday! 

65  5  1 

66  4  ! 

67  4  2 

68  5  1 

69  4  2 

70  5  3 

71  6  3 

72  J*  2 

73  3  2 

7^  ^  2 

75  4  1 

76  J*  2 

77  ^  2 

78  3  1 

79  ^  2 

80  3  1 

81  2  1 

82  3  1 

83  2  1 

84  1  0 

85  1  1 

Total  75  32~ 


Table  6.6  -  Weekday  and  Weekend  Temperature  Distribution  for 
Model  Summer. 
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Teraperature  (°F)      Weekdays      Weekends/Holidays 

65  4  2 

64  2  2 

63  2  1 

62  2  1 

61  1  0 

60  2  1 

59  2  1 

58  2  0 

57   2  1 

56  2  0 

55  1  o 

5^  2  1 

53  2  l 

52  l  l 

51  2  1 

50  2  2 

Total  31  nr 


Table  6.7  -  Weekday  and  Weekend  Temperature  Distribution  for 
Model  Fall. 
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weekend/holiday  respectively.   Referenced  to  the  period 
January  1976  to  February  1977,  the  steam  consumption  for  any 
day  is  a  function  only  of  outside  ambient  temperature. 

Assignments  of  daily  total  kilowatt  demand  were  made  in 
accordance  with  the  percentage  distributions  derived  in 
Chapter  IV.   For  the  representative  year  the  number  of  week- 
days and  weekends  falling  within  each  of  the  temperature 
ranges  depicted  in  Figures  4.3  -  4.11  was  determined.   This 
number  was  multiplied  by  each  of  the  bar  graph  percentages 
which  show  the  relative  proportion  of  days  in  each  tempera- 
ture grouping  that  fall  within  specified  five-kilowatt 
bandwidths.   (The  lower  value  of  kilowatt  demand  was  used  in 
each  bandwidth.)   The  above  procedure  resulted  in  a  propor- 
tionate dispersion  of  daily  total  kilowatt  loads  for  days 
within  specific  temperature  ranges  of  the  model  year.   Actual 
association  of  a  daily  total  demand  with  one  particular  tem- 
perature day  was  made  randomly.   Although  this  method  lacks 
the  sophistication  which  characterized  the  steam  load  assign- 
ment, it  ensures  a  normative  spread  of  electrical  consumptions, 
consistent  with  1976  and  early  1977  data. 

With  both  steam  and  electrical  24  hour  demands  thus 
enumerated  for  the  model  year,  daily  load  profile  assignments 
were  made. 

6. 5   Daily  Profile  Assignment 

The  steam  demand  profiles  which  were  developed  in 
Chapter  III  for  each  season  were  applied  directly  to  weekdays 
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and  weekends/holidays.   It  is  recalled  thai-  the  demand  model 
included  two  load  patterns  for  winter  and  spring  weekdays, 
one  being  an  "extreme"  profile.   Assignment  of  this  profile 
to  particular  days  within  the  winter  and  spring  seasons  was 
made  randomly,  the  only  restriction  being  that  the  overall 
percentage  of  weekdays  showing  this  pattern  be  30%. 

Electrical  load  profiles  were  prescribed  according  to 
the  frequency  of  occurrence  in  the  sample  year  (see  page  99  )  . 
Individual  assignments  were  made  randomly  so  as  to  achieve 
the  correct  proportionate  distribution  of  profiles  for  each 
temperature  band. 

For  purposes  of  computer  simulation  of  MIT  demands,  a 
numbering  sequence  was  used  to  denote  each  profile  for 
efficient  decision-making  use  within  the  simulation  program. 
The  tvoe  of  dav  was  specified  as  the  variable  TYPDAY.   Steam 
profiles  were  denoted  bv  the  variable  STMPRO.   Electrical 
profiles  were  described  bv  KWPRO.   A  listing  of  the  variables 
and  the  numerical  designations  for  each  appears  below: 

Variable    Designation  within  Program   Meaning 

TYPDAY  Type  of  Day  Being 

Modeled 

1  Weekday 

2  Weekend/Holidav 

STMPRO  Steam  Profile  Number 

1  Normal  Winter 

2  Extreme  Winter  or 

Spring 

3  Normal  Spring 

4  Normal  Summer 

5  Normal  Fall 
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variable   Designation  within  Program   Meaning 

KWPRO  Electrical  Profile 

Number 

1  Weekday  with  Extreme 

Mid-Afternoon  Peak 
(T>60°F) 

2  Weekday  with  Normal 

Mid-Afternoon  Peak 
(T>60°F) 

3  Weekday  with  Morning 

Peak  (T>60°F) 

4  Weekday  with  Late 

Afternoon  Peak 

(T  £  60°F) 

5  Weekday  with  Early 

Afternoon  Peak 

(T  £  60°F) 

6  Weekend  with  3:00 

P.M.  Peak 

7  Weekend  with  5:00 

P.M.  Peak 

8  Weekend  with  6:00 

P.M.  Peak 


6  .  6   I n t egration  of_ Model  Year  Data  Into  a  Computer 
Program 

A  method  is  desired  which  permits  the  simple  transfer 

of  load  information  into  a  simulation  program.   The  description 

which  follows  is  purposefully  general  in  that  the  modeling  of 

different  total  energy  system  configurations  may  require 

slight  modifications  of  the  main  program  sequence. 

6.6.1   Data  Input 

In  the  simulation  program,  each  steam  and 
electrical  profile  is  input  as  a  single  array  of  twenty-four 
elements.   Within  each  array  the  respective  hourly  load  factors 
are  listed  sequentially  from  1:00  A.M.  to  midnight.   The 
designations  for  the  arrays  are  as  follows: 
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Steam  Profiles 

WWD  Winter  weekday 

XWD  Extreme  Winter  or  Spring  Weekday 

WWE  Winter  Weekend 

SPWD  Spring  Weekday 

SPWE  Spring  Weekend 

SUWD  Summer  Weekday 

SUWE  Summer  Weekend 

FAWD  Fall  Weekday 

FAWE  Fall  Weekend 

Electrical  Profiles 

WDA1  Weekday  above  60°F,  Profile  #  1 

WDA2  Weekday  above  60  °F,  Profile  #  2 

WDA3  Weekday  above  60°F,  Profile  #  3 

WDB1  Weekday  below  60°F,  Profile  #  1 

WDB2  Weekday  below  60°F,  Profile  #  2 

WEN1  Weekend  Profile  #  1  (3:00  P.M.  Peak) 

WEN2  Weekend  Profile  #  2  (5:00  P.M.  Peak) 

WEN3  Weekend  Profile  #  3  (6:00  P.M.  Peak) 

Each  day  of  the  model  year  is  input  separately,  fully 
described  by  a  listing  of  five  numbers  on  a  data  card.   The 
numbers  are  the  values  assigned  to  the  following  constants 
within  the  program: 

TEMP  Average  ambient  temperature  for  day  (°F) 

TYPDAY  Weekend  or  Weekday 

STMPRO  Steam  profile  number  assignment 

KWPRO  Electrical  profile  number  assignment 

DAYKW  Daily  total  kilowatt  load  (as  explained  in 

Section  6.4) 

6.6.2   Program  Sequence 

The  program  is  designed  to  read  all  the  pro- 
file arrays  first.   After  initializing  several  parameters  for 
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later  use  in  the  simulation,  it  then  reads  the  first  day's  data. 
Based  on  the  information  in  TYPDAY  and  STMPRO,  the  program 
stores  the  proper  steam  load  profile  in  the  array  HSLF  (Hourly 
Steam  Load  Factor)  .   Similarly,  the  number  assigned  to  KWPRO 
governs  the  assignment  of  the  correct  electrical  profile  to 
HELF  (Hourly  Electrical  Load  Factor) . 

Equations  3.1  and  3.2  are  included  as  statement  functions 
within  the  main  program.   Dependent  upon  whether  the  day  being 
simulated  is  a  weekend  or  weekday,  a  computation  is  made  of 
daily  total  steam  demand  (DAYSTM) .   Daily  kilowatt  demand  is 
read  as  input  data  (DAYKW) . 

Individual  subroutines  may  be  designed  to  model  specific 
total  energy  system  configurations.   Any  number  of  arguments 
may  be  specified  in  the  designation  of  each  subroutine,  but  as 
a  minimum  the  following  four  must  be  included:   HELF,  HSLF, 
DAYKW,  DAYSTM.   This  ensures  that  the  necessary  load  informa- 
tion for  one  day  is  transferred  to  each  subroutine  where  it 
will  be  simulated  as  hourly  steam  and  electrical  demands  to 
specific  pieces  of  equipment. 

It  is  envisioned  that  each  subroutine  will  contain  pro- 
visions for  the  calculation  of  fuel  consumption  and  waste 
heat  available  and  then  return  this  information  to  the  main 
program.   The  main  program,  in  turn,  may  be  designed  to  keep 
a  running  record  of  the  fuel  consumed,  kilowatt  demand,  etc. 
In  addition  it  may  be  structured  to  compute  the  appropriate 
cost  of  power  purchased  from  Cambridge  Electric  Company  for 
the  case  when  the  particular  total  energy  system  is  not 
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supplying  all  of  MIT's  electrical  needs  (see  Chapters  VII  and 
VIII).   Thus,  much  flexibility  in  modeling  is  possible. 

The  program  continues  by  reading  a  second  day's  data. 
It  repeats  the  above  procedures  for  passing  information  to 
each  subroutine.   As  more  days  are  input  with  different  steam 
and  electrical  demands,  certain  total  energy  system  designs 
will  begin  to  appear  more  attractive  than  others. 

Figure  6.1  is  a  generalized  flow  chart  for  the  manipula- 
tion of  input  data  within  the  program. 

6.6.3   Arrangement  of  Data  Deck 

The  365  days  of  input  data  are  grouped  by 
season.   They  are  not,  however,  arranged  sequentially  by 
increasing  (or  decreasing)  temperature  within  any  one  season. 
Rather,  an  attempt  has  been  made  to  group  days  near  the  mean 
historical  monthly  temperatures  for  the  Boston  area  (see 
Table  6.8).   The  intent  is  to  gain  as  much  realism  in  the 
model  as  possible  so  monthly  cost  breakdowns  will  be  meaning- 
ful. 

The  fact  that  the  representative  year  included  a  greater 
number  of  warmer  and  colder  days  than  the  historical  smoothed 
temperature  distribution  prevented  a  direct  grouping  about 
the  monthly  mean  temperatures  of  Table  6.8.   Instead,  tempera- 
tures were  biased  lower  in  winter  and  higher  in  summer  to 
account  for  the  expanded  temperature  distribution  of  Table  6.3, 
Typically,  days  within  10 °F  of  the  monthly  average  were  chosen 
as  being  representative  of  the  monthly  temperature  spread. 
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Month  of  Year 

Average 

Normal  Temperature 

January 

29.2 

February 

30.4 

March 

38.1 

April 

48.6 

May 

58.6 

June 

68.0 

July 

73-3 

August 

71.3 

September 

64.5 

October 

55.^ 

November 
December 

45.2 
33.0 

Table  6.8  -  Historical  Monthly  Average  Temperature 
for  Boston  (1941  -  1970). 
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Table  6.9  shows  the  final  monthly  temperature  averages  used  in 
the  data  deck  arrangement. 

The  value  of  the  above  procedure  is  quickly  recognized 
if  an  attempt  is  made  to  interpret  operating  costs  over 
periods  of  time  which  are  smaller  than  a  year.   One  energy 
option  which  is  open  to  MIT,  perhaps  unlikely,  is  the 
installation  of  a  partial  electrical  generation  capability 
with  the  balance  of  electrical  needs  (peaking)  furnished  by 
Cambridge  Electric.   Billing  by  Cambridge  Electric  is  accom- 
plished on  a  30  day  basis.   If  the  grouping  of  days  within 
any  one  season  is  not  made  according  to  some  average  monthly 
temperature,  the  months  of  June,  July  and  August  could 
possibly  show  very  similar  electrical  bills.   This  would  be 
misleading  and  detracts  from  the  credibility  of  the  load  model. 

A  listing  of  representative  year  temperature  and  load 
information  has  been  included  as  Table  6.10.  Data  is  shown 
sequentially  by  month  for  the  365  days  of  the  model  year.   It 
may  be  used  directly  in  a  simulation  program. 

6 . 7  Need  for  Validation  of  Representative  Year  Model 
If  total  energy  system  operating  costs  are  to  be 
compared  with  existing  costs  for  providing  campus  energy,  the 
modeling  of  specific  plant  configurations  must  be  accomplished 
using  load  information  typical  of  some  "representative"  year 
at  MIT.   This  has  been  the  objective  of  the  work  thus  far. 
The  detailing  of  daily  steam  and  electrical  usage  patterns 
so  as  to  reflect  typical  operating  conditions  is  complete.   A 
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Month  of  Year  Average  Temperature  (°F) 
January  26.7 
February  28.0 
March  35-7 
April  ^8.6 
May  58.6 
June  71 . 1 
July  76.5 
August  74.^ 


September 


November 


December 


66.^ 


October  55- ^ 

1*2.8 


30.6 


Table  6.9  -  Average  Daily  Temperature  for  Months 
in  Representative  Year  Model. 
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means  of  ascertaining  the  validity  of  the  load  model,  however, 
is  required.   To  this  end,  Chapter  VII  has  been  devoted. 
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VII   VALIDATION  OF  LOAD  MODEL 

Prior  to  use  of  the  demand  model  in  evaluating  alternative 
total  energv  svstem  designs,  it  is  necessarv  to  establish  its 
validitv.   One  method  of  determininq  the  accuracy  of  load 
estimation  is  to  compare  the  loads  predicted  by  the  model  with 
those  actually  observed  during  some  previous  time  period  for 
which  reliable  temperature  information  is  available.   As  tem- 
perature is  the  only  parameter  of  interest,  it,  alone,  may  be 
used  as  a  yard-stick  for  measuring  the  worth  of  the  load 
model. 

7. 1   Simulation  of  Central  Utility  Plant  Operation 

Temperature  data  for  the  calendar  year  1976  was  used 
in  a  computer  program  designed  to  simulate  the  operation  of 
MIT's  existing  Central  Utility  Plant.   The  program  was  con- 
structed according  to  the  generalized  flow  chart  of  Chapter  VI. 
Specified  for  each  day  of  the  year  were  outside  average 
temperature  (TEMP) ,  daily  steam  and  electrical  load  profile 
designations  £TMPRO  and  KWPRO) ,  and  daily  total  kilowatt 
demand  (DAYKW) .   A  subroutine  called  STMMIT  was  written  for 
the  purpose  of  simulating  the  fuel  consumption  characteristics 
of  the  200  psig  boilers  at  the  Central  Utility  Plant.   Figure 
7.1  is  a  listing  of  the  program  and  subroutine.   A  description 
of  program  variables  follows: 

Main  Program 

TOTKW  -  a  storage  location  which  starts  at  0.0  for 
each  month  and  which  sums  the  daily  total 
kilowatt  load  for  a  30  day  period. 
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TOTFUE  -  storage  location  which  sums  the  daily  fuel 
consumed  for  a  30  day  period. 

M  -  a.  counter  which  indexes  by  one  (1)  for  each 
day  being  simulated. 

K  -  a  counter  which  indexes  up  to  30,  reflecting 
load  simulation  for  one  month,  used  as  a 
criterion  for  printing  output. 

J  -  a  counter  which  corresponds  to  the  month 
for  which  loads  are  beinq  simulated. 

HR170  -  designation  for  long  hours  use;  special 
charges  apply  to  kilowatt  usage  above  a 
set  level  for  each  month. 

HRLONG  -  the  dollar  charge  for  long  hour  use  by 
Cambridge  Electric. 

ENCHRG  -  energy  charge  for  purchased  electricity. 

DMDCHG  -  demand  charge  for  purchased  electricity. 

RATEAD  -  utility  rate  adjustment  (presently  15.6%) 

FUELAD  -  fuel  adjustment  which  accounts  for  increased 
cost  of  fuel  purchased  by  utility  (presently 
2.549C  per  KWh) . 

COSTKW  -  total  monthly  cost  of  purchased  electricity. 

Subroutine  STMMIT 

AVGSTM  -  average  hourly  steam  demand. 

AVGKW  -  average  hourly  electrical  demand. 
HRELEC  -  hourly  electrical  demand. 

HRSTM  -  hourly  steam  demand. 

PEAKHR  -  peak  electrical  demand  for  24  hours;  main 
program  has  provision  for  storing  peak 
monthly  demand. 

DMD1  -  hourly  demand  for  boiler  #3  at  Central  Plant. 
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DMD2  -  hourly  demand  for  boiler  #4  at  Central  Plant. 
DMD3  -  hourly  demand  for  boiler  #5  at  Central  Plant. 

FUEL1  -  fuel  consumed  in  satisfying  DMD1. 

FUEL2  -  fuel  consumed  in  satisfying  DMD2. 

FUEL3  -  fuel  consumed  in  satisfying  DMD3 . 
FUEL  -  daily  total  of  fuel  consumed. 
STMMIT  causes  hour  by  hour  steam  loads  to  be  imposed  on 
mathematical  models  of  boilers.   In  reality,  no  more  than 
three  boilers  are  ever  used  (under  present  campus  loads)  to 
satisfy  demand  requirements.   Therefore,  instructions  were 
provided  in  STMMIT  to  simulate  one  boiler  operation  for  loads 
less  than  70,000  lbs  steam/hour,  two  boiler  operation  for 
loads  less  than  140,000  lbs  steam/hour  and  three  boiler 
opeation  for  all  loads  greater  than  this.   Boiler  fuel  consump- 
tion rates  were  obtained  from  Mr.  George  Reid  at  the  Central 
Utility  Plant  and  incorporated  into  program  statements.   With 
each  hourly  steam  demand  imposed  on  the  plant,  the  subroutine 
calculated  the  amount  of  fuel  oil  (#6  residual)  required  to 
satisfy  the  demand.   For  each  day  being  simulated  a  24  hour 
total  of  fuel  consumed  was  returned  to  the  main  program. 

A  provision  for  monitoring  the  daily  peak  electrical  load 
was  incorporated  into  STMMIT.   The  purpose  here  was  to  simu- 
late monthly  billing  by  Cambridge  Electric.   Daily  kilowatt 
demands  for  the  model  year  (DAYKW)  were  used  in  place  of  1976 
load  information  to  determine  the  predicted  electrical  costs 
for  a  typical  year  at  MIT  under  the  present  rate  structure. 
Demand  charges,  energy  charges  and  long  hour  charges  were 
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computed  in  the  main  program  for  each  30  day  period.  A  rate 
adjustment  of  15.6%  and  a  fuel  adjustment  of  2.549C/KWh  were 
applied  in  accordance  with  current  billing  regulations  under 
Rate-8. 

7.2   Steam  Load  Model  Results 

The  results  of  the  simulation  were  most  encouraging. 
As  a  predictor  of  steam  demand,  the  model  proved  excellent. 
Figure  7.2  displays  1976  consumption  information  relative  to 
that  estimated  by  the  program  based  on  1976  temperatures. 
It  can  be  seen  that  only  one  model  month  deviated  substantially 
from  1976  consumption  levels.   The  program  overpredicted  the 
September  steam  demand  by  approximately  18%.   Two  explanations 
were  found  for  this  anomoly. 

During  the  summer  months,  two  chiller  units  are  typically 
used  to  satisfy  campus  cooling  requirements.   It  has  been  the 
recent  practice  of  the  operating  personnel  at  the  Central 
Utilitv  Plant  to  cease  operation  of  the  larger  of  these  units 
durinq  the  fall  season  to  conserve  fuel.   Since  coolinq  demand 
is  lower  at  this  time  than  durinq  the  summer,  one  unit  is 
usuallv  sufficient  to  handle  the  requisite  loads.   It  does 
happen,  however,  that  davs  occur  in  the  fall  for  which  the 
cooling  demand  cannot  be  satisfied  by  one  chiller  unit  alone. 
Nevertheless,  a  second  unit  might  not  be  placed  on  line  for 
the  simple  reason  that  such  swings  in  outside  ambient  tempera- 
ture are  transitory.   It  is  not  a  prudent  engineering  practice 
to  constantly  shift  chiller  units  on  and  off  for  the  sake  of 
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ensuring  that  the  mix  of  equipment  is  optimum  for  the  parti- 
cular load.   Although  an  obvious  answer  lies  in  running  two 
chiller  units  during  the  fall,  each  at  a  reduced  load  factor, 
such  operation  is  wasteful  from  an  efficiency  standpoint. 
There  is  no  set  policy  which  prohibits  two  chiller  units  from 
remaining  in  operation  into  the  fall;  depending  upon  the 
daily  temperature  trend,  this  course  of  action  might  possibly 
be  followed.   It  happened  in  1976,  however,  that  September 
was  a  moderate  month  and  two  chiller  operation  was  avoided 
as  much  as  practicable.   The  result  was  a  reduction  in  the 
amount  of  steam  used  by  the  Chiller  Plant  relative  to  that 
which  would  have  been  consumed  under  two  unit  operation. 

A  second  reason  for  lower  overall  fuel  consumption  during 
September  concerns  the  reduction  in  steam  supplied  to  portions 
of  the  main  qroup  for  heatinq  purposes.   Under  normal  circum- 
stances, a  20  inch  header  conducts  5  psiq  steam  from  the 
exhaust  of  each  central  plant  turbine  driven  auxiliary  to 
areas  of  the  main  group.   During  the  summer  months  this 
header  is  closed  off.   No  firm  date  in  the  fall  exists  for 
its  reopening.   Rather,  this  decision  depends  upon  the  neces- 
sity for  heating  in  campus  buildings.   Typically,  motor 
driven  auxiliary  equipment  is  used  during  warm  weather  opera- 
tion of  the  plant  as  there  is  no  use  for  the  low  pressure 
steam  which  exhausts  from  turbine  driven  auxiliaries.   As  the 
need  for  heating  arises  (into  the  fall  season) ,  the  20  inch 
header  is  opened  and  a  gradual  shift  to  turbine  driven  equip- 
ment is  made.   Although  there  were  days  during  September  of 
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1976  for  which  heating  would  have  been  desirable,  opening  of 
the  header  was  delayed  until  decidedly  "winter"  weather  could 
be  foreseen.   As  a  result,  a  tendency  to  operate  only  the 
electric  auxiliaries  prevailed.   Overall  steam  consumption 
was,  therefore,  lower  than  it  might  otherwise  have  been. 

For  the  entire  year  the  computer  model  overpredicted 
fuel  consumption  by  1.4%.   Six  of  the  months  showed  less  than 
.1%  difference  between  that  which  was  actually  consumed  and 
that  which  was  estimated.   The  greatest  disparities  were  noted 
in  the  spring  and  fall  seasons,  presumably  because  of  peculiari- 
ties associated  with  the  shift  from  heating  to  air  conditioning 
and  back  again. 

Temperature  data  for  the  representative  year  (Chapter  VI) 
was  used  in  a  second  computer  run.   As  model  year  temperatures 
vary  from  those  in  1976,  it  was  desired  to  determine  how  well 
monthly  fuel  consumption  figures  reflect  the  temperature 
difference.   Figure  7.3  summarizes  information  relative  to  fuel 
consumption  and  temperature  in  1976  and  the  model  year.   The 
temperature  plots  are  expressed  relative  to  the  historical 
monthly  averages  for  Boston.   It  is  verified  that  consumption 
levels  track  temperature  closely.   Figure  7.4  is  the  computer 
output  for  the  simulation  of  representative  year  demands.   It 
should  be  noted  that  fuel  totals  were  computed  on  the  basis 
of  a  30  day  month.   Prior  to  inclusion  in  Figures  7.2  and  7.3 
they  were  scaled  to  the  appropriate  number  of  days  in  each 
month. 
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7 . 3   Electrical  Load  Model  Results 

Comparison  of  1976  and  representative  year  electri- 
cal demands  revealed  that  the  temperature  model  accounts  for 
the  major  differences  between  the  two  sets  of  monthly  consump- 
tion figures.   Relative  to  1976  kilowatt  totals,  months  in 
the  representative  year  showed  lower  winter  and  higher  summer 
electrical  demands.   Figure  7.5  illustrates  the  data  trend. 
This  difference  is  a  result  of  the  distribution  of  temperature 
in  the  model  year.   It  may  be  recalled  from  Section  6.2.1  that 
refinement  of  the  historical  temperature  model  resulted  in  an 
expansion  of  both  upper  and  lower  temperature  extremes.   For 
the  representative  year  model,  assignment  of  daily  total 
kilowatt  loads  was  made  in  accordance  with  the  percentage 
distributions  of  Chapter  IV.   In  that  the  monthly  average 
temperatures  for  the  model  year  deviate  substantially  in  some 
cases  from  1976  averages  (see  Figure   7.3),  it  is  to  be 
expected  that  kilowatt  demand  should  also  differ. 

Figure  7.4  shows  that  the  peak  electrical  demand  predicted 
by  the  load  model  is  15,583  KW.   While  this  is  higher  than  the 
most  recent  peak  of  15,240  KW,  it  is  not  significant  enough  to 
warrant  revision  of  the  entire  model.   Inclusion  of  this  peak 
will  ensure  a  conservative  proposed  total  energy  plant  sizing. 
The  fact  that  a  purely  random  application  of  the  various 
daily  electrical  load  profiles  resulted  in  a  peak  so  close 
to  the  present  high  is  indicative  of  the  quality  of  the 
electrical  load  model. 
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The  1976  yearly  total  electricity  purchase  from  Cambridge 
Electric  was  85,296,000  KWh.   As  predicted  by  the  representa- 
tive year  temperature  assignment,  the  annual  MIT  electrical 
consumption  should  be  83,600,000  KWh.   It  is  clear  that 
although  the  load  model  apportions  the  consumption  differently 
by  month  (in  accordance  with  temperature  distribution) ,  the 
net  result  in  terms  of  total  usage  is  very  similar.   It  can 
be  said,  therefore,  that  1976  was  not   atypical  from  a  stand- 
point of  demand. 

The  average  cost  of  purchased  electricity,  as  calculated 
in  the  computer  program  for  the  model  year  (Figure  7.4), 
agrees  almost  perfectly  with  current  cost  figures  available 
through  the  offices  of  the  Physical  Plant.   Numbers  on  the 
order  of  3.64£/KWh  are  typical  of  recent  billings  from 
Cambridge  Electric.  Until  such  time  as  a  new  rate  adjustment 
or  fuel  adjustment  is  authorized,  the  sequence  of  program 
steps  concerned  with  computation  of  monthly  billing  charges 
will  remain  valid. 

7 . 4   Further  Application  of  Demand  Model 

On  the  basis  of  the  preceding,  simulation  of  specific 
total  energy  schemes  may  be  undertaken.   Campus  growth  can  be 
allowed  for  as  a  simple  percentage  addition  to  the  representa- 
tive year  demands.   That  is,  daily  total  kilowatt  and  steam 
loads  may  be  multiplied  by  a  constant  to  simulate  any  magni- 
tude of  demand  increase.   The  validity  of  the  load  model  has 
been  established.   It  remains  to  devise  a  methodology  for 
modeling  particular  total  energy  designs. 


-195- 

VIII  SELECTION  OF  PLANT  DESIGN  &  METHODOLOGY  FOR  MODELING 

The  possible  choices  of  equipment  configuration  for  the 
proposed  total  energy  system  design  are  numerous.   The  most 
practical  schemes,  however,  center  on  three  general  types  of 
plants;  steam  extraction,  gas  turbine  and  diesel.   Within  any 
particular  plant  classification,  an  abundance  of  design 
variations  exist.   The  addition  of  helper  turbines  and  waste 
heat  boilers  to  individual  cycles,  for  example,  affords  a  high 
degree  of  flexibility  to  certain  plant  designs.   A  detailed 
examination  of  the  many  possible  engineering  alternatives  is 
not  intended.   Rather,  an  overview  of  the  three  general  design 
arrangements  is  envisioned  with  particular  attention  devoted 
to  outlining  the  methods  of  modeling  each  for  computer 
simulation. 

8 . 1   Steam  Extraction  System 

The  use  of  a  steam  extraction  system  at  MIT  would 
require  installation  of  higher  pressure  boilers  than  those 
which  presently  exist  at  the  Central  Plant.   For  the  size  of 
generation  facility  in  question  (21  MW) ,  pressures  on  the 
order  of  800  psig  are  typical.   It  is  likely  that  condensing 
turbines  which  have  a  single  automatic  extraction  capability 
would  be  best  suited  for  MIT's  needs.   An  extraction  pressure 
of  200  psig  would  provide  the  requisite  process  steam  for 
campus  heating  and  the  Central  Chiller  Plant  use. 

The  sequence  of  operation  for  a  single  automatic  steam 
extraction  system  is  straightforward.   High  pressure  steam  is 
provided  to  turbines  which  drive  electrical  generators.   Steam 
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is  exhausted  to  a  condenser  at  a  vacuum  of  approximately 
3"  Hg  A.   A  variable  amount  of  steam  may  be  extracted  from  the 
turbine  at  a  constant  pressure,  independent  of  the  flow  to  the 
turbine  driven  electrical  generator.   In  the  event  that  the 
flow  of  extraction  steam  is  not  sufficient  to  fulfill  demand 
requirements,  the  steam  supply  may  be  augmented  through  a 
reducing  valve  off  the  high  pressure  steam  main  or  from  the 
existing  low  pressure  boilers.   Boiler  firing  rate  is  a  func- 
tion of  the  composite  demands  for  electrical  power  generation 
and  extraction  steam. 

For  MIT's  purposes  a  single  generator  sized  to  accommo- 
date all  electrical  needs  is  feasible  although  it  does  not 
afford  a  backup  capability.   The  use  of  two  generators,  each 
appropriately  sized,  together  with  the  newer  two  or  three 
existing  boilers,  offers  the  advantage  of  flexibility  in 
operation.   For  the  sake  of  illustrating  how  an  extraction 
system  may  be  modeled,  single  generator  operation  has  been 
assumed  for  MIT. 

8.1.1   Mathematical  Model 

Information  from  the  manufacturer  on  steam 
turbine  generator  performance  is  typically  in  the  form  of 
straight  line  graphs.   Figure  8.1  is  an  example  of  the  perfor- 
mance characteristics  for  a  General  Electric  15,000  KW 
generator  with  a  single  automatic  extraction  at  200  psig. 
For  any  particular  electrical  demand  the  graph  indicates  what 
range  of  steam  extraction  is  available.   For  a  chosen 
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extraction  level,  throttle  flow  (boiler  demand)  may  be  read 
on  the  y-axis.   It  is  observed  that  while  the  generator  is 
sized  for  15,000  KW,  loads  of  25,000  KW  are  achievable  (thus, 
appropriately  sized  for  MIT) .   With  increasing  load  factor, 
the  minimum  extraction  flow  also  increases.   The  possibility 
exists,  therefore,  that  more  extraction  steam  might  be  avail- 
able than  is  needed  for  heating/air  conditioning  purposes. 

In  order  to  translate  the  information  contained  in 
Figure  8.1  into  a  series  of  program  statements  designed  to 
model  the  system  operation,  a  schematic  is  required.   Prior 
to  outlining  this  procedure,  however,  the  equations  of 
several  lines  on  the  graph  are  needed.   The  following  descrip- 
tion is  purposefully  general  so  that  any  particular  steam 
extraction  turbine  may  be  modeled,  providing  only  that  perfor- 
mance lines  similar  to  those  in  Figure  8.1  are  available. 

The  equation  of  the  line  describing  extraction  as  a 
function  of  load  for  the  maximum  throttle  flow  must  be 
determined.   An  easy  method  of  accomplishing  this  is: 

(a)  Extend  the  line  of  maximum  throttle  flow  toward  the 
right  in  Figure  8.1. 

(b)  Extend  several  of  the  lines  of  constant  extraction 
until  they  intercept  the  line  drawn  in  (a) . 

(c)  Read  the  magnitude  of  generator  output  at  each 
intersection  of  the  lines  from  parts  (a)  and  (b) . 

(d)  Plot  a  straight  line  graph  of  extraction  versus 
generator  output  using  the  information  from  (c) . 
Determine  the  slope  (nO  and  y-intercept  (b^  of 
the  line: 
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ymax   xtr   =  mlx  +  bl  (Equation    I) 

where 

x  =  generator  output  (15  MW  a   x  £  25.5  MW) 

y  =  maximum  possible  extraction  at  specified 
generator  output 

Figure  8.2  illustrates  the  above  procedure.   It  should  be 
noted  that  Equation  I  is  valid  only  for  generator  loads  falling 
within  15  and  2  5.5  MW. 

A  second  equation  which  is  needed  is  that  for  minimum 
permissible  extraction  as  a  function  of  generator  output, 
applicable  within  the  same  range  of  loads  as  indicated  above 
(15  -  25.5  MW) .   By  reading  the  magnitudes  of  generator  output 
at  each  intersection  of  the  maximum  exhaust  flow  line  and  the 
lines  of  constant  extraction,  a  straight  line  plot  of  generator 
load  versus  minimum  extraction  may  be  drawn.   Figure  8.3 
illustrates  the  procedure.   The  slope  (nu)  and  y-intercept 
(b2)  of  the  line  may  be  easily  determined  with  the  result  that 
an  equation  of  the  following  form  is  constructed: 

y  .       =  nux  +  b„  (Equation  II) 

■^min  xtr     2     2  ^ 

where 

x  =  generator  output  (15  MW  <  x  <  25.5  MW) 

y  =  minimum  extraction  at  specified  generator  output 

For  generator  loads  less  than  the  rated  capacity  (15  MW) 
an  expression  relating  maximum  permissible  extraction  to 
generator  output  is  required.   The  intersection  of  constant 
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Intersec- 
tion  of 
maximum 
exhaust 
flow  line 
and  lines 
of  constant 
extraction. 

(2)  Read  magni- 
tudes below 
and  plot 
line. 


Figure  8.3 


Generator  Output 

Illustration  of  Method  for  Determining  the 
Equation  of  Minimum  Permissible  Extraction 
as  a  Function  of  Generator  Load  (Valid  for 
Loads  Greater  Than  the  Rated  Generator 
Capacity) . 


(1)  Intersection  of  constant 
extraction  lines  with 
line  of  minimum 
exhaust  flow. 

(2)  Read  magnitudes 
below  and  plot 
line  of  maximum 
extraction. 


Generator  Output 

Pigure  8.4  -  Illustration  of  Method  for  Determining 
the  Equation  of  Maximum  Permissible  Ex- 
traction as  a  Function  of  Generator  Load 
(Valid  for  Loads  Less  Than  Rated  Capacity) 
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extraction  lines  with  the  line  of  minimum  exhaust  flow  defines 
the  points  needed  to  construct  an  appropriate  straight  line 
plot.   Valid  only  for  loads  less  than  15  MW,  Equation  III  is 
determined  in  the  same  manner  as  were  the  previous  two  with 
slope  (m  )  and  y-intercept  (b_)  easily  computed: 


where 


ymax  ext  =  m3X  =  b3  (Equation  III) 


x  =  generator  output,  <  15  MW 

y  =  maximum  permissible  extraction  at  specified 
generator  output 


Figure  8.4  highlights  the  above  procedure. 

The  last  equation  which  must  be  derived  for  the  purpose 
of  modeling  is  that  which  yields  throttle  flow  as  a  function 
of  generator  load.   It  may  be  expressed  as 


m.x  +  b.  (Equation  IV) 


where 


y  =  throttle  flow 
m.  =  slope  of  lines  of  constant  extraction 

x  =  generator  output 
b.  =  b'  +  minimum  extraction  intercept 

The  minimum  extraction  intercept  is  defined  as  the  extension 
of  the  zero  extraction  line  to  where  it  intersects  the  zero 
generator  output  ordinate.   It  may  be  visualized  as  the 
hypothetical  minimum  throttle  flow  for  an  unloaded  generator, 
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although  this  interpretation  is  for  explanatory  purposes  only. 
The  value  of  b1  depends  upon  the  amount  of  steam  extracted  at 
any  particular  time.   It  is  defined  by  the  following  relation: 

b»     XTR 


Ab     AXTR 


where 


Ab  =  maximum  extraction  intercept  -  minimum  extraction 
intercept 

XTR  =  amount  of  steam  being  extracted 

AXTR  =  maximum  permissible  extraction  -  minimum 
permissible  extraction 

=  maximum  H.P.  turbine  flow  -  minimum  L.P. 
turbine  flow 

The  maximum  extraction  intercept  may  be  viewed  as  the  hypo- 
thetical maximum  throttle  flow  for  an  unloaded  generator,   b* 
represents  a  proportionate  increase  in  throttle  flow,  over  the 
hypothetical  minimum  at  zero  extraction  which  results  from  an 
increase  in  the  demand  for  extraction  steam.   Consideration 
of  it  arises  because  of  the  different  scales  depicted  on 
Figure  8.1  for  extraction  and  throttle  flow.   Figure  8.5 
illustrates  the  above  description. 

8.1.2   Program  Schematic 

The  following  sequence  of  steps  permits 
organization  of  a  computer  program  designed  to  simulate  the 
operation  of  a  single  automatic  extraction  steam  turbine 
generator.   It  has  been  assumed  that  both  steam  and  electrical 
loads  have  been  placed  on  the  system,  and  it  is  required  to 
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determine  the  requisite  boiler  steam  flow  which  will  satisfy 
campus  demands. 

(a)  Determine  if  the  electrical  load  falls  within  the 
rating  of  the  specific  generator  (0  -  25.5  MW) .   If 
the  load  is  greater  than  the  upper  limit,  instruc- 
tions must  be  provided  to  divide  the  demand  between 
two  separate  units  (using,  perhaps,  a  small  diesel 
generator  to  provide  excess  generating  capacity) . 

(b)  Assuming  the  electrical  demand  falls  within  the 
correct  range,  ascertain  whether  it  is  greater  or 
less  than  the  generator  rating  (15  MW) .   If  it  is 
less,  skip  to  step  (h) . 

(c)  If  the  electrical  load  is  greater  than  the  generator 
rating,  use  Equation  I  to  determine  whether  the 
steam  demand  can  be  supplied  totally  by  extraction 
steam.   (If  the  steam  demand  is  greater  than  the 
maximum  permissible  extraction  at  that  generator 
load,  steam  supply  must  be  augmented  by  reducing 
high  pressure  steam  or  by  utilizing  the  existing 
boilers. ) 

(d)  For  steam  demands  which  are  less  than  the  maximum 
permissible  extraction,  it  must  be  verified  that  all 
the  steam  which  is  extracted  can,  in  fact,  be  used. 
Equation  II  is,  therefore,  applied  to  find  the 
minimum  permissible  extraction. 

(e)  If  the  minimum  allowable  extraction  is  greater  than 
the  steam  demand,  some  of  the  excess  steam  must  be 
dumped  as  wasted  heat.   The  alternative  would  be 

to  operate  the  generator  at  a  load  which  just 
satisfies  the  required  extraction  steam  demand  and 
purchase  the  balance  of  electricity  from  Cambridge 
Electric.   Thus,  a  feature  should  be  included  in  the 
program  to  track  the  quantity  and  frequency  of  such 
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mismatches.  It  is  suggested  that  the  most  efficient 
means  of  accomplishing  this  would  be  through  the  use 
of  a  separate  subprogram. 

(f)  If  the  steam  demand  proves  to  be  greater  than  the 
minimum  extraction,  it  is  verified  that  all  the 
steam  which  is  extracted  can  be  used.   Throttle  flow 
is,  therefore,  computed  in  accordance  with  Equation 
IV  with  XTR  equal  to  the  steam  demand. 

(g)  In  the  event  that  it  is  not  necessary  to  supplement 
the  flow  of  extraction  steam  to  satisfy  demand 
(Step  c) ,  boiler  load  is  equal  to  throttle  flow. 
Otherwise,  it  is  equal  to  throttle  flow,  computed 
with  XTR  equal  to  maximum  extraction,  plus  the 
amount  of  required  augmenting  steam. 

(h)   For  the  case  when  the  electrical  demand  is  less  than 
the  design  generator  rating,  a  determination  must  be 
made  as  to  whether  the  steam  demand  is  greater  or 
less  than  the  maximum  extraction  at  that  electrical 
load.   Equation  III  is  entered  with  a  value  for 
required  generator  output.   It  yields  the  value  of 
maximum  permissible  extraction  flow.   A  simple 
comparison  between  the  magnitude  of  this  number  and 
that  of  steam  demand  indicates  whether  augmenting 
steam  is  necessary. 

(i)   Throttle  flow  (boiler  load)  is  computed  in  accordance 
with  Equation  IV  with  XTR  equal  to  steam  demand, 
assuming  the  demand  does  not  exceed  the  maximum 
allowable  extraction.   Otherwise,  XTR  equals  maximum 
extraction  in  Equation  IV,  and  boiler  demand  becomes 
the  sum  of  throttle  flow  and  augmenting  steam. 

For  any  variation  of  a  steam  extraction  total  energy 
icheme  the  above  instructions  may  be  followed.   Use  of  this 
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sequence  of  steps  in  the  construction  of  a  simulation  program 
ensures  an  accurate  performance  model  for  a  particular  genera- 
tor sizing.   By  combining  boiler  performance  information  with 
the  boiler  demands  which  result  from  the  above  modeling 
procedure,  an  appreciation  for  the  relative  cost  advantages 
of  a  steam  extraction  total  energy  system  can  be  gained. 

8.  2   Diesel  Generator  System 

Consideration  of  a  diesel  generator  system  for 
supplying  MIT's  electrical  needs  implies  that  at  least  some  of 
the  existing  boilers  at  the  Central  Utility  Plant  would  remain 
in  operation.   The  alternative  would  be  the  complete  use  of 
unfired  boilers  to  recover  some  of  the  heat  rejected  from  the 
diesel  engines.   For  MIT's  purposes,  it  is  envisioned  that  a 
diesel  generator  system  with  waste  heat  boilers  could  be 
economically  feasible.   Although  the  design  pressure  of  such 
auxiliary  units  might  not  be  high  enough  to  supply  normal 
campus  steam  demands,  feed  heating  operations  could  be 
supported  by  the  waste  heat  boilers.   They  would  serve  as  a 
supplement  to  the  newer  existing  boilers  (#3,  4,  and  5), 
therby  increasing  plant  efficiency.   A  model  can  be  developed 
to  simulate  the  hourly  consumption  of  waste-heat  steam,  pro- 
viding a  method  for  measuring  plant  efficiency  increases.   It 
is  suggested  for  preliminary  design  estimates,  however,  that 
since  the  existing  boilers  alone  will  supply  the  campus  steam 
demand,  modeling  of  the  overall  power  plant  may  be  simplified 
by  not  considering  the  waste  heat  units  directly.   Their 
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influence  may  be  taken  into  account  as  a  percentage  reduction 
in  fuel  consumption  over  a  range  of  Central  Utility  Plant 
boiler  loads.   The  discussion  which  follows  assumes  the  use 
of  the  present  MIT  boilers  in  combination  with  several  diesel 
generators. 

The  number  of  diesel  generators  which  could  conceivably 
be  incorporated  into  a  total  energy  plan  is  a  function  of 
availability  of  space,  rated  generator  capacity  and  the 
requirements  for  backup  power.   While  a  cost  tradeoff  study 
will  likely  narrow  down  the  available  options,  the  modeling 
of  a  diesel  configuration  can  proceed  in  the  absence  of 
definitive  generator  sizes  and  numbers.   Once  a  simulation 
scheme  is  devised,  system  parameters  may  be  easily  altered, 
providing  only  that  appropriate  performance  curves  are  avail- 
able for  each  generator  sizing. 

Performance  information  for  a  diesel  generator  is  typi- 
cally provided  from  the  manufacturer  as  a  curve  of  brake 
specific  fuel  consumption  versus  load.   Unlike  that  for  the 
steam  extraction  turbine,  the  diesel  qenerator  performance 
information  is  sliqhtlv  nonlinear.   Construction  of  a  math 
model,  therefore,  follows  a  different  approach  than  was  used 
in  the  previous  section. 

As  a  first  step  in  the  procedure,  the  equation  for  the 
performance  curve  of  each  unit  must  be  determined.   Least 
square  regression  techniques,  like  those  described  in  Chapters 
III  and  IV,  may  be  utilized  for  this  purpose.   Similarly, 
linear  interpolation  procedures  are  effective,  providing  a 
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sufficient  number  of  reference  points  are  used.   Once  a. 
method  is  defined  which  permits  the  determination  of  fuel 
consumption  as  a  function  of  generator  load,  the  following 
procedures  are  recommended. 

(a)  Model  steam  demands  in  a  manner  analagous  to  that 
used  for  the  plant  simulation  in  Chapter  VII.   As 
the  calculation  of  boiler  fuel  consumption  is 
entirely  independent  of  electrical  load  fluctuation, 
a  separate  subroutine  may  be  used  to  model  the 
boiler  operation. 

(b)  Specify  as  a  variable  the  number  and  rated  capacities 
of  the  diesel  generators.   A  series  of  instructions 
should  be  provided  for  the  (simulated)  start  up  of 

a  second  generator  once  the  load  on  the  first 
reaches,  for  example,  80%  of  its  rated  capacity. 
An  efficient  way  of  accomplishing  this  is  through 
program  statements  similar  to  those  in  the  STMMIT 
subroutine  of  Chapter  VII  (Figure  7.1).   Depending 
upon  the  rated  capacity  of  the  generators  being 
modeled,  simulated  start  up  or  shut  down  of  an  unit 
could  occur  every  few  hours  as  loads  fluctuate  during 
the  day.   Some  care,  therefore,  must  go  into  specify- 
ing generator  operating  limits  to  avoid  an 
unrealistic  simulation  model. 

(c)  For  each  diesel  in  operation  a  calculation  of  hourly 
fuel  consumption  should  be  made  as  loads  vary  during 
the  day.   By  incorporating  the  equation  of  each 
performance  curve  into  a  statement  function  within 
the  program,  electrical  demands  may  be  conveniently 
satisfied  and  appropriate  fuel  costs  computed. 

The  actual  operation  of  a  diesel  generator  system  at  MIT 
would  most  likely  deviate  from  the  above  model  in  the  manner 
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in  which  alternate  units  are  placed  on  line.   As  opposed  to 
automatic  start  up  of  a  second  generator  at  a  preset  load,  the 
central  computer  of  the  Facilities  Management  System  would 
generate  an  advisory  message  to  plant  personnel  based  on  the 
predicted  campus  demand.   Depending  upon  the  rate  of  load 
increase,  start  up  of  a  second  generator  could  occur  in 
advance  of,  for  example,  a  preset  80%  load  on  the  first.   For 
determining  the  relative  cost  differences  of  alternate  total 
energy  designs,  however,  the  procedure  outlined  yields  valid 
consumption  information.   It  is  not  anticipated  that  overall 
plant  efficiency  would  be  markedly  affected  by  a  model  which 
ignores  the  potential  management  benefits  of  the  FMS  computer. 

8 . 3   Gas  Turbine  Configuration 

Installation  of  a  system  of  gas  turbine  driven 
electrical  generators  at  MIT  would  most  likely  be  accompanied 
by  the  addition  of  waste  heat  boilers  to  the  Central  Plant. 
Depending  upon  the  capacity  of  each  waste  heat  unit  and  the 
provisions  for  firing  it  separately,  one  or  more  of  the 
existing  Central  Plant  boilers  might  be  needed  for  supplemen- 
tary generation  of  steam.   For  the  size  of  gas  turbines 
available  for  use  at  MIT,  however,  waste  heat  boilers  with 
sufficient  generation  capacity  to  satisfy  peak  winter  demands 
are  manufactured  commercially.   Whether  MIT  would  elect  to 
dispose  of  its  present  complement  of  boilers  in  order  to  pro- 
vide the  necessary  space  for  several  waste  heat  boiler 
additions  is  a  question  which  can  only  be  resolved  after  a 
design  study  is  completed.   For  discussion  sake  a  methodology 
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for  modeling  an  integral  gas  turbine/waste  heat  boiler  total 
energy  design  is  addressed. 

Like  the  diesel  generator,  gas  turbine  performance 
characteristics  are  summarized  by  a  single  curve.   When  the 
system  includes  an  exhaust  boiler,  the  manufacturer  provides 
information  on  steam  flow  by  way  of  a  second  curve  which  is 
defined  over  the  range  of  gas  turbine  loads.   Additionally, 
a  separate  boiler  performance  curve  exists  for  use  when  the 
unit  must  be  auxiliary  fired.   Prior  to  constructing  a  system 
model,  an  analytic  expression  for  each  performance  curve  is 
required.   As  gas  turbines  have  characteristically  poor  fuel 
consumption  at  off  design  loads,  their  performance  curves  are 
non  linear.   Consequently,  linear  regression  or  linear  inter- 
polation procedures  are  recommended  for  the  purpose  of 
devising  a  means  of  representing  the  performance  information 
in  a  programmable  manner.   A  suggested  scheme  for  system 
modeling  follows. 

(a)  Input  as  program  variables  the  number  and  maximum 
rated  capacity  of  each  gas  turbine/waste  heat  boiler 
configuration. 

(b)  Develop  a  series  of  program  instructions  which 
assign  particular  units  to  be  in  operation  for  a 
range  of  specified  electrical  loads.   Satisfaction 
of  steam  requirements  will  be  a  secondary  considera- 
tion. 

(c)  For  each  level  of  electrical  load  placed  on  the 
generator (s)  compute  the  fuel  required  to  power 
the  gas  turbine  and  the  amount  of  steam  which  can 
automatically  be  furnished  at  that  load  from  the 
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respective  waste  heat  boiler (s) . 

(d)  Compare  the  steam  demand  with  that  available  from  the 
waste  heat  boilers.   If  the  demand  exceeds  what  is 
being  provided,  the  boilers  must  be  separately  fired. 
The  program  should  be  designed  so  as  to  keep  a 
record  of  the  frequency  of  need  for  additional 
steam.   For  each  such  instance,  a  computation  of  the 
fuel  consumed  must  be  made. 

(e)  If  campus  steam  demand  is  less  than  the  supply  from 
the  unfired  boilers,  the  program  should  store  infor- 
mation on  the  degree  of  mismatch.   Although  the 
steam  generation  rate  may,  in  reality,  be  reduced  by 
venting  some  of  the  gas  turbine  exhaust,  output 
should  be  available  to  the  program  user  on  the  per- 
centage of  the  time  campus  demand  was  lower  than 
the  design  level  of  steam  generation. 

As  was  the  case  with  the  diesel  generator  modeling 
procedure,  actual  operation  of  a  gas  turbine/waste  heat  boiler 
system  differs  from  the  description  above.   Operating  personnel 
must  anticipate  the  need  for  augmenting  the  generation  of 
steam.   The  model,  however,  assumes  no  time  lag  from  when  the 
demand  is  perceived  until  it  is  satisfied.   By  incorporating 
the  features  of  FMS  into  the  system  operation,  as  would  likely 
occur  at  MIT,  switching  of  gas  turbines,  for  example,  in 
advance  of  a  predictable  peak  could  occur.   Nonetheless,  a 
computer  model  which  affords  a  considerable  amount  of  flexi- 
bility in  plant  simulation  may  be  constructed  in  the  absence 
of  FMS  considerations. 
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8. 4   Overview  of  Modeling  Procedure;   Cost  Analysis 

While  a  comparison  of  alternate  total  energy  system 
designs  can  be  made  on  the  basis  of  relative  differences  in 
annual  fuel  consumption  (plant  efficiency) ,  several  other 
factors  exert  strong  influence  upon  a  design  feasibility 
study.   Features  may  be  incorporated  into  a  computer  program 
to  account  for  these  additional  modeling  considerations.   The 
end  result  is  a  better  appreciation  of  the  relative  cost 
advantages  of  one  design  over  another. 

The  most  crucial  yardstick  for  plant  comparison  is  that 
of  acquisition  cost.   This  includes  not  only  the  purchase  of 
machinery  but  also  the  costs  related  to  installation  and  test- 
ing of  equipment.   Additionally,  it  covers  the  initial  monetary 
outlay  for  spare  parts.   In  the  face  of  growing  pressure  to 
install  pollution  abatement  devices  on  all  power  plants,  MIT 
will  witness  sizable  acquisition  cost  increases  for  any 
proposed  total  energy  system. 

Annual  operating  costs  are  an  important  basis  for  com- 
parison of  different  system  designs.   Over  and  above  the 
expenditures  for  fuel,  maintenance  costs  are  included  in  this 
category.   Plant  insurance  and  manning  costs  are  recurring 
annual  expenses  as  well.   While  plant  acquisition  cost 
estimates  proceed  largely  on  the  basis  of  quotes  from  the 
manufacturer,  operating  costs  are  more  difficult  to  predict. 
The  uncertainty  in  fuel  prices  provides  a  sizable  threshold 
for  error  in  a  comparative  plant  study. 

A  convenient  means  of  placing  the  forementioned  costs  in 
the  proper  perspective  is  through  a  model  of  monetary 
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expenditures  based  on  the  annual  cost  method  of  accounting. 
It  is  assumed  that  the  capital  required  for  procurement  of  a 
total  energy  system  would  come  through  the  sale  of  public 
bonds  by  MIT.   The  following  equations  describe  the  time 
stream  of  payments  for  an  annual  cost  evaluation. 
Annual  monetary  outlay  by  MIT: 

M  =  C   x  (R/P  )  +  C  ($/yr) 

a       o     o 


where 


and 


/   _   i  x  (l+i) 

7°~    d+i,n-i 


R  =  total  annual  payment  covering  the  interest  on 
the  sale  of  bonds. 

i  =  interest  rate  of  bonds. 

P   =  amount  of  money  received  as  a  result  of  the  bond 
sale  (lump  sum) . 

n  =  number  of  years  for  amortization. 

C   =  acquisition  cost, 
a 

C   =  yearly  operating  cost. 

It  is  envisioned  that  a  20  year  payment  plan  would  be  chosen 
as  this  represents  the  approximate  life  of  a  new  power  plant. 

Some  means  of  depreciating  the  chosen  total  energy  system 
must  be  included  in  a  cost  model.   Since  MIT  is  a  non-profit 
making  organization,  consideration  must  be  given  to  the 
annual  savings  which  result  from  supplying  on-site  electrical 
power  as  opposed  to  purchase  from  Cambridge  Electric.   For  the 
purpose  of  computer  modeling,  it  should  be  assumed,  therefore, 
that  a  certain  proportion  of  the  annual  savings  is  set  aside 
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each  year  to  account  for  depreciation.   With  this  money 
conservatively  invested,  enough  should  be  available  after  20 
years  to  support  purchase  of  new  equipment.   While  this  is  not 
what  actually  occurs  for  capital  investments  at  MIT,  inclusion 
of  a  program  feature  to  account  for  depreciation  provides  a 
further  measure  of  the  economic  attractiveness  of  a  chosen 
total  energy  design. 

The  cost  analysis  should  center  on  a  comparison  between 
the  equivalent  cost  of  purchased  electricity  from  Cambridge 
Electric  and  the  annual  costs  for  sustaining  a  total  energy 
system  in  operation.   The  yearly  depreciation  should  be  sub- 
tracted from  the  equivalent  cost  of  purchased  electricity. 
Positive  savings  result  when  the  annual  costs  to  MIT  for  pro- 
viding its  own  electricity  are  less  than  the  equivalent 
purchased  electrical  costs  (minus  the  depreciation) . 

8.5   Summary 

Analysis  of  plant  performance  is  fundamental  to  a 
comparative  study  of  alternative  means  for  providing  MIT's 
energy  needs.   Several  total  energy  system  designs  have  been 
outlined  for  possible  use  at  MIT.   A  methodology  for  carrying 
out  a  numerical  simulation  of  the  operating  conditions  for 
each  has  been  presented.   By  varying  the  specific  mix  of 
equipment,  through  statement  changes  in  a  computer  program, 
a  wealth  of  information  can  be  gained  on  the  likely  choice  of 
a  total  energy  system  for  MIT. 
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